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PREFACE 

The procedures and methods presented herein are intended primarily to 
serve a3 a guide and text to those whose work consists of developing smooth, 
continuous lines with a uniformly changing rate of curvature. These requirements 
are met in many engineering fields of endeavor, including the development of 
aircraft and automobile shapes, ship lines, projectile and missile shapes and 
trajectories, highway and railroad curves, and many other similar branches of 
engineering. 

The increased speeds required of airplanes, ships, missiles, railroads, 
etc., - to mention only a few cases - necessitate extreme care being paid to 
the design of their external shapes so that proven drag penalties, both in 
operating cost and performance, may be minimized. Hence it is no longer suffic- 
-ient, or advisable, to resort to the old method of drawing these lines "by eye,* 
but instead these lines must be developed from analytical expressions which 
will insure that uniform rates of change in curvature are obtained, and that 
even the minutest reversals in curvature are eliminated. Both of the above condit- 
-ions must be satisfied in order to keep the drag of any given shape as small as 
possible. 

In the further case of railroad track and highway layout, the smooth 
transition from a straight track, or highway, into a curve and then back again 
into a straight course, requires that these transitions be started and ended 
gradually, and that the rate of change of curvature vary smoothly along the entire 
length. Here again analytical expressions defining these curves are needed to 
insure that these requirements are satisfied. 

The use of such analytical expressions for the shapes of many bodies 
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also permits the easy calculation of the volumes, surface areas, etc., of these 

I 

bodies, and allows proportional changes in shape, volume, or area, to be quickly 

made. Thi3 reduces considerably the time required making successive layouts, 
each time checking to determine if the desired conditions have been obtained. 

The development of equations for such purposes as outlined above, and for 
many other similar cases which will suggest themselves to the reader, is given in 
detail herein. The equations have been kept as simple as possible, with only 
elementary algebra and calculus employed, so that it is not necessary for the 
reader to be well versed in either of these subjects in order to make use of the 
equations, or develop those of his own. The text, however, will serve as a handy 
guide and reference to elementary calculus, as it represents a current practical 
application of the calculus. 

These methods of analytically developing curved lines and streamlined 
shapes have been employed for many years in the design of aircraft manufactured 
by the Republic Aviation Corporation, whose airplanes have long been recognized 
for their exceptionally high performance. No small amount of this performance 
gain is due to the extreme care taken in, the analytical development of the shapes 
involved, so that smooth airflow may be maintained, and drag losses kept to a 
minimum. 

In this regard, the author wishes to express his deep appreciation of 
the efforts of Mr. A. Kartveli, Chief Engineer of Republic Aviation Corporation, 
whose insistence upon perfection in design led to the need for the development 
of the methods and equations given on the following pages. Grateful acknowledg- 
-ment is also extended to Dr. W. J. O'Donnell and Messrs. Hoshen Lu and J. Cravero, 
of Republic Aviation Corporation, for their aid and recommendations in the 
preparation of the manuscript. 
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FOREWORD 

The generation of faired, curved surfaces is a problem frequently 
encountered in such engineering applications as ship building, airplane construct- 
-ion, automobile body design, etc. This generation is usually performed graphically, 
and is based on the geometrical. properties of known curves, usually of the second 
degree. In many cases, surfaces are even generated without any reference to known 
curves simply by "feel and eye*. 

In the course of the last few years, however, the rigid requirements 
of aeronautical design made it necessary to devise a much more accurate analytical 
method for the generation of curved surfaces. This method consists in performing 
the work of generation purely analytically by representing various curves evolved 
by algebraic equations. By a combination of such equations, desired shapes can be 
obtained with a great degree of accuracy. The use of equations, in addition, enables 
the computation of areas, volumes, centers of gravity, etc., to be considerably 
simplified. 

Mr. Harry Haase has been engaged in this type of work for several years 
at the Republic Aviation Corporation, and thiy book represents the results of his 
experience. His method of expressing the ordinates of a curve as a function of the 
abscissae, by means of a polynomial with fractional exponents, proves to be very 
successful in generating wings and streamlined bodies, and in fairing the intersect- 
-ion between two bodies. In fact, the efficiency of the method is so great that it 
has now obsoleted and replaced other methods previously used by the Engineering 
Department of the Republic Aviation Corporation. 

This book will bring a valuable contribution to the Engineering liter- 
-ature on this subject, and should prove a very useful guide to designers engaged 
in this type of work. 

A. Kartveli 

Chief Engineer 

Republic Aviation Corporation 



CHAPTER I 
THE CIRCLE AND THE ELLIPSE 

1.00 - TTsr of the Ctrnle in the Development n f Contour Lines:- In the 
development of contour lines for bodies moving through a fluid, such as 
air or water, two types of lines are of prime interest to the designer. 
These are the 'horizontal plane contour lines, or lines parallel to the 
direction of fluid flow which are usually called streamlines, and the lines 
perpendicular to the direction of fluid flow, or cross-sectional lines. 
Contour lines in other planes are generally not necessary when the above- 
mentioned contour lines are thoroughly and accurately developed. 

The circle is of particular importance in the development of 
streamlines , or cross-sectional contour lines. While it is seldom, if ever, 
used directly as a streamline, it is often used in the development of 
certain airplane wing airfoil sections, and in laying out the leading edges 
of streamlined shapes. As a cross section, the circle is roost efficient for 
a body moving through a fluid, and is expensively used for airplane fuselages, 
rockets, missiles, pressurized vessels, etc. 

A shape with circular cross sections along its entire length is 
known as a "body of revolution", which is formed by rotating any line, curved 
or straight, 360 degrees about its- reference axis. Hence, rotating a straight 
line completely about a reference axis parallel to it generates a cylinder, 
the simplest type of body of revolution, while rotating, for example, the 
upper contour line of an airfoil about its reference axis will produce a 
streamlined body of revolution. 

If a volume of fluid is forced past a body of revolution parallel 
to its reference axis, or if the body is pushed or pulled through a stationary 
fluid, the fluid displaced by this body will move, at a given cross section 
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or station, at the same speed, and the pressure will be -uniform around the 
periphery for each station. Therefore, no cross-flow can result at any 
station. The fluid disturbances thus created will be small, and the drag 
produced will then be a minimum. 

In the case of any other cross section, such as an elliptical one, 
for example, the fluid flowing along the major axis will travel farther than 
that moving along the minor axis, since the distance is longer due to the 
increased curvature. The fluid will then increase in speed along the minor 
axis to keep pace with the free fluid, but will, increase in speed more along 
the major axis for the same reason. A difference in pressure will then result 
around the periphery at a given station, and a cross-flow created to cause 
fluid disturbances and consequently more drag. The ellipse, however, has its 
place in the development of various contour lines, as will be shown later. 
1.10 -Equations of the Circler- 

^ l.H — Origin at Center ;- The Pythagorean theorem states that, in a 
right triangle, the square of the length of the hypotenuse is equal to the 
sum of the squares of the lengths of the two other sides. By use of this 
theorem, it is possible to locate any point on the periphery of a circle of 
a fixed radius if either the abscissa or the ordinate is known, or, if the 
radius is not known, it can be determined when the abscissa and ordinate of 
a point are given. 

In the circle shown in Fig. 1, let (p) be any point on the periphery 

x 




FIGURE 1 
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at a distance (x) from the vertical reference line, or Y axis, and distance 
(7) from the horizontal reference line, or X axis, and let (r) be the radius. 
Then, from the above theorem: 

y 2 +x 2 =r 2 (1) 

. or: 

y . \/r2 _ x2 ...(2) 

The slope at any point on the periphery may be found by means of the 

first derivative of Eq. (1): -:;•:-. 

2y(dy/dx) ■»- 2x = 

(dy/dx) - -(x/y) . ....(3) 

The second derivative of Eq. (1), which is often used in line 
development work for the purpose of determining points of inflection, or 
change in curvature, is found by differentiating the first derivative as given 
in Eq. (3): 

(d2y/dx 2 ) = (Vy 2 )(dy/<ix) .--- 1/7 "-.;-_ .. 

- (x/y2)(-x/y) - 1/V 

- -(*2/y3 + 1/y) 

.= -(i/y 5 )(x 2 + y 2 ) -". — ( 4 > 

Substituting from Eq. (1): 

(d2y/dx2) = -(r2/y3) .(4a) 

Very often it is easier to use the trigonometric equations of a 
circle in layout work, such as when laying out circles of large radius. In 
Fig. 1, for example, it is evident that: «, 

(y/r) m sin 6 

y = r(sin e) (5) 

and; 

(x/r) s cos 9 

x = r(cos 9) ( 6 ) 
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From Eqs. (5) and (6) the ordinates of a circle of any given radius 
(r) can be easily plotted by assuming values of and computing, with the 
aid of trigonometric tables, the corresponding values of (x) and (y). 

1.12 - Transfer of Origin ;- In the previous equations the origin of 
the circle was assumed to be at the center. It is more convenient, though, 
in some instances to have the origin on the X axis and on the periphery, as 
at point (B) in Fig. 2. This is particularly true when another curve is added 
to an arc of a circle, as in the case of a leading edge of an airfoil section. 
All abscissae can then be measured from one point, the transferred origin of 
this circle. 




FIGURE 2 
Since the distance from the center of the circle, (0), to point (B) 
is equal to the radius (r), the Pythagorean theorem of Eq. (1) can be written! 

y2 + (x - r)2 - r 2 
or: 

y 2 + x 2 - 2rx + r2 = r2 



2rx - x2 



(7) 



y » ferx - x2 (7a) 

The first derivative, or slope at any point on the periphery, is: . 

(dy/dx) = < r - x > .-..(8) 

\JZrx - x^ 

Substituting y =\/2rx - x 2 from Eq. (7a) in Eq. (8), the slope is: 
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(dy/dx) = (r - x)/y (9) 

The second derivative is obtained as follows by differentiating 
Eq. (9): 

(dVdx 2 ) = [-D(r - xhT^dy/dx) - y* 1 
Rearranging, and substituting the value of (dy/dx) from Eq. (9): 

(dVd* 2 ) - - ((r ■- x)/y 2 )((r - x)/y) - (l/y) 
= - (r-x)2/y3_ ( 7 2/ 7 3) 

." ' ■ - (-(r - x) 2 -7 2 )/7 3 ..............(10) 

Eqs. (7a), (9), and (10) only apply when the origin is placed on the 
X axis at the periphery of the circle, as shown in Fig. 2. In this form, the 
development of faired lines can be most easily accomplished, and the equations 
are simpler to use. If the origin of the circle should be displaced both from 
the X axis and the periphery, as in Fig. 3, the general equation of the circle 
becomes: 
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FIGURE 3 









(x - h) 2 + (y _ k) 2 = r 2 (11) 

The derivatives of Eq. (11) become considerably more complicated to 
use than those of Eq. (7a), and since it is seldom, if ever, necessary in line 
development work to transfer the origin in such a manner as shown in Fig. 3, 
Eq. (11) need not be further discussed. 
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Transfer of the origin to a point not on the periphery, but on the 

X axis, is, however, useful when it is necessary to add to the left hand 

side of the circle. This is shown in Fig. 4 below: 

T 




FIGURE 4 



In this case, it is readily seen that the equation of the circle 
is expressed as: 

yS+Cx-fa^-r 2 ....(12) 

or: 

^ :■;:?■>"■■ y 2 + x 2 - 2hx + h 2 > r 2 ■,"■":.. 

and solving the above for (y): 

~~ (13) 



y ss\|r 2 - x 2 - h 2 +■ 2hx 

1.13 - ?r,lnr (inordinates of the Circle ;- The equations of a circle can 
also, be expressed in polar coordinate form, instead of the rectangular 
coordinate fo~n previously discussed. This is shown in Fig. 5, where the 
origin, or pole, is assumed to be on the periphery of the circle and also on 
the X axis. 




FIGURE 5 
Again using the theorem that the square of the hypotenuse of a right 
triangle is equal to the sum of the squares of the other two sides: 

^ = x2+ 7 2 -- .........(14) 

and, from trigonometry: 

x = ^>)cos(e) ...........(15) 

The basic expression for transfer of the origin is given by Eq. (7), 

which may also be written: 

x 2 + y 2 - Erx (16) 

Hence, if (/>) 2 is substituted in Eq. (16) for (x 2 + y 2 ), and cos © 
for (x), the polar equation of a circle with the origin on the X axis and on 
the periphery becomes: 

r> 2 = 2r( cos 9) 

or: J3 = 2r(cos 6) (17) 

For polar coordinates of a circle with the origin any distance (h) from 
the center, as shown in Fig. 6, but still on the X axis, Eqs. (14) and (15) 
still apply. 

Substituting the values of (x 2 + y 2 ) and (x) from Eqs. (14) and (15), 
respectively, in the general equation of a circle with the origin transferred 
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FIGURE 6 
any distance (h) along the X axis, (Eq. (12a)), we obtain the final equation: 

jC?2 = r 2 - h 2 + 2h^>cos 9 (18) 

1.20 - aflwpt.nt.1on nf t.hP TCI n pap t.n DpvpI nnmRnt. of Contour Linen:- The ellipse 
is one of the roost widely used geometric figures in the development of contour 
lines, being extensively used for cross-sections and segments of streamlined 

shapes. 

As a cross-section of an airplane fuselage, for example, it has not 
only been used for convenience in designing the encloseure of a cockpit, leav- 
ing ample and suitable space in the fuselage for fuel tanks, air ducts, controls, 
equipment, and the like, but it has a definite advantage over a circular cross- 
section in the case of a mid-wing airplane, where the wing panels join the 
fuselage. It was stated before that the fluid flowing along the major axis of 
an ellipse has to travel a longer distance than the fluid moving along the 
minor axis. It will, therefore, have to move faster along the top of the 
fuselage, ( it is assumed that the major axis of the ellipse is vertical), to 
keep up with that flowing along the sides. 

However, the wing too has thickness, and the fluid travelling along 
the minor axes of the ellipses of the fuselage will be forced to go over or 
under the wing surfaces with increased velocity. The fluid speeds will thus be 
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equalized, to some extent, in the region up to the maximum thickness of the 
airfoil section. From there on the air will expand, and the ratio of major 
to minor axes of the ellipses should be gradually reduced until they finally 
end up in a circle at the rear of the fuselage. If the horizontal tail 
surfaces are attached directly to the fuselage, an elliptical cross-section 
should be maintained, since the same velocity changes will occur. 

In aircraft work, another item of interest may be noted with regard 
to the location of the wings with respect to the fuselage. The upper surface 
of conventional airfoil sections usually has more curvature than the lower 
surface, and therefore a higher fluid velocity occurs over the upper surface. 
This velocity differential, when expressed by means of Bernoulli^ equation 
in terms of its corresponding pressure differential, is necessary for the 
creation of lift. To more nearly equalize the fluid velocities over the wing 
and fuselage, therefore, it is advisable to locate the wing slightly below the ■ o 
centerline of the fuselage. This helps to preserve a smooth flow at the 
junction, and has been found by experience to materially reduce the so-called 
"interference" drag. Obviously, the amount of displacement of the wing from 
the fuselage centerline should be determined from a consideration of whether 
the fuselage cross-section is circular, elliptical, or otherwise, and also on 
a consideration of the airfoil section being used. 
1.30 - Equations and Construction of the Ellipse t- 

1.31 _ Standard Equation of the Ellipse :- The ellipse, which in reality 
is a flattened circle with uniformly varying curvature, can be constructed by 
the use of two circles with a common center at the origin. The radius of the 
larger circle becomes the major axis of the ellipse, while the radius of the 
smaller circle is the minor axis. 

Vectors passing through the origin may then be drawn at random. For 
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each vector a line can be drawn parallel to the X axis passing through the 
intersection of this vector with the periphery of the small circle. Similarly, 
for this same vector, a line can be drawn parallel to the Y axis passing 
through the intersection of this vector with the periphery of the larger 
circle. The intersection of these two lines then determines one point on the 
periphery of the ellipse, such as point (p) on Fig. 7. Continuing this 
process, as many points as desired can be found, and joined together to form 
the ellipse. 




FIGURE 7 
It is evident from Fig. 7 that: 

sin 8 s (x/a) •••' 

and: 

cos 9 = (y/b) * ' 

Using the relationship (sin 2 + cos 2 ) = 1, from trigonometry, Eqs. 
(19) and (20) may be combined to form the standard equation of an ellipse: 

(x/a) 2 + (y/b) 2 = 1 



(19) 



( 2 0) 
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or: 



x 2 /a 2 + 7 2 /b 2 = 1 ... 
Eq. (21) can be solved for (y) as follows: 
y 2 = b 2 ((l - (x 2 /a 2 )) 
= (b 2 /a 2 )(a 2 - x 2 ) 
or: y - (b/a) \/a 2 - x 2 ... 



.(21) 



.(22) 



The slope of the curve at any point on the periphery is determined, 
as in the case of the circle, by finding the first derivative of the equation 
for the curve. Hence, using Eq. (21) and differentiating: 

2x/a 2 + 2y(dy)/b 2 (dx) = 
from whence: 

(dy)/(dx) - -OAOAaV) < 23 > 

The second derivative, which is used for determining the rate of 

change of slope, or the locations of any points of inflection, is found as 

follows: 

(dy)/(dx) - -(bM/CaV) * 
= -(b 2 x)(a 2 y)~ 1 

Then: 

(d^/dx 2 ) = (-b 2 x)(-l)(a 2 y)(a 2 )((dy)/(dx)) + (a^)' 1 ^ 2 ) 

=((b 2 xa 2 )/(a 2 y) 2 )(dy/dx) - (b 2 )/(aV) 

>,^ 2 (dv/dx) - b 2 (a 2 y) 
(a^) 2 

Substituting the value of (dy/dx) from Eq. (23): 

= -b 2 (b 2 x 2 + a 2 y 2 )/(a 4 y 5 ) 

If Eq. (21) be multiplied by (a 2 b 2 ), there results (b^ 2 ) 4- (a^ ) 
(a 2 b 2 ). Therefore, substituting this equation in the above: 



12. 

(d*y)/(dx 2 ) = -((b 2 (aV))/(aV 3 ) 

- -(b^/CaV 5 ) (24) 

Eq. (24) then gives the value of the second derivative of the 

standard equation of an ellipse. 

1.32 - Plotting an Ellipse bv Means of Trigo nometric Functions:- A 
convenient method of plotting an ellipse is by means of trigonometric tables. 
By assuming any value of (x), and letting (a) be the major axis, the sine of 
the angle can be computed by employing Eq. (19). The cosine can then be 
interpolated from trigonometric tables, and (y) found by multiplying the 
cosine by (b), the minor axis. This establishes a point on the periphery. Repeat- 
ing this process for other values of (x), a number of points can be plotted, 
and connected by means of a spline or french curves. The angles are not needed 
in degrees, and do not have to be noted. 

An example of the method of computing the ordinates of an ellipse 
with a major axis of (a) - 10 and minor axis of (b) - 6 is shown in Table I: 

TABLE I 



X 


(x/a) a sin 6 


cos 9 


y = (b) cos 


2 


0.200 


0.979 


5.874 


4 


0.400 


0.916 


5.496 


5 


0.500 


0.866 


5.196 



Plotting the results, a curve as shown in Fig. 8 will be obtained: 



X 




Assumed y's I 
i 



*■ Assumed x's 
FIGURE 8 
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It will be found that It becomes quite difficult to plot an accurate 
ellipse when the angle 6 reaches its maximum. Therefore, the procesB of 
computing the ordinates can be reversed by assuming values of (y) and then 
calculating the corresponding values of (x). It is advisable, for a good check, 
to compute the (y's) up to about 40 degrees and then start- with values of (x) 
from about 30 degrees so that the points overlap. The Table will be the same 
as Table I except that the first column will be (y), the second (y/b) » cos e » 
the third sin 6, and the fourth (x) - (a) sin 9. It is obvious that only one 
quadrant need be computed, as the others can be copied. 

1.35 _ Transfer of Origin ;- Just as in the case of a circle, it Is 
often more convenient to transfer the origin, to a point on the X axis and on 
the periphery of the ellipse. 

Transferring the origin from (0) to (B), as shown in Fig. 9, the 

distance between these two points is equal to (a), the major axis of the 

T 
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FIGURE 9 
ellipse. Employing the standard equation of the ellipse, as given by Eq. (21), 
and substituting (x - a). for (x), we get: 

(x - a)Va 2 + y 2 A 2 - 1 
(x 2 - 2ax + a 2 )/a2 + y2/b 2 = 1 
Multiplying both sides of the equation by a 2 b 2 : 
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or: 



b 2 (x 2 - 2ax + a 2 ) + a 2 ^ = a 2 b 2 

y 2 * (2ab 2 x - b^J/a 2 
- (b 2 /a 2 )(2a 3 c-x 2 ) 



and: 



_ (b 2 /a 2 ) [" (a 2 b 2 - Sab^ 4- h 2 ^) - a 2 ^ 2 ] 
L a 2 ? -3 J 



-a 2 b 4 + Sab^x - b 4 x 2 - a^ 2 ? 2 
aV 5 



(27) 



y > (b/a) teax - x 2 (25) f 

It will be noted that Eq. (25) becomes the same as Eq. (7a), the i 

' I 
equation of a circle with its origin transferred as in this case, when I 

• •'•■■ ■ ' * ; fi 

b - a = x, the radius of a circle. I 

The slope at any point on the periphery, (dy/dx), is found as follows: f 

" I 

From Eq. (25): | 

y - (b/a)(2ax - x 2 )^ | 

Then: 1' 

(dy/dx) = (b/a)(J)(2ax - x 2 H(2a - 2x) 
« ((b(a - x))/ a \/2ax - x 2 
Substituting the value of V2ax - x 2 from Eq. (25): 
(dy/dx) - b(a - x)/(a)(ay/b) 

- b 2 (a - xj/a^- ..(26) 

The second derivative in this case is: 
(dSy)Acbc 2 ) = (b 2 /a 2 ) [(a - x)(-y- 2 (dy)/(dx» + (y" 1 - 1)] 

- (b*/a 2 ) £-((a - x)(dy/dx))/y 2 - 1/y] 
Substituting the value of (dy/dx) from Eq. (26): 
(d^/dx 2 , -(b 2 /a*) [(' (a 1 x)(b 2 /a 2 ) ((a - x)/r))/y 2 + 1/y} 
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By letting a - b = r, the radius of a circle, in Eq. (27), the 
expression for the second derivative of the equation of a circle with its 
origin transferred to the periphery on the X axis, Eq. (10), is obtained, as 
would be expected. 

1.34 - Polar Coordinates of an Ellipse ;- The. polar coordinates of an 
ellipse are sometimes used in layout work, although perhaps not as often as 
the rectangular coordinates due to slightly more required computation. Once 
the necessary calculations have been made, however, the actual layout is some- 
what simpler as it is only needed to mark off various distances on a given set 
of radii. 

In Fig. 10, below, we can see that: 

I 




T 
FIGURE 10 



x = R(sin ^) 



and: 



y » R(cos <f) 
where (R) is the radius to any given point (P) on the periphery of 

the ellipse. 

If these values be substituted in Eq. (21), the standard equation of 

an ellipse in rectangular coordinates, then: 
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(R(sia *0) 2 /a 2 - (R(cos jrf)) 2 /b 2 = 1 
or: 

R 2 (sinVa 2 - cosVb 2 ) = l 
This can be reduced to: 

R - (ab)/ l/b 2 sinV - a 2 cos 2 ^ ( 28) 

Eq. (28) then gives an expression for the radius to any point on the 
periphery of the ellipse in terms of the major and minor axes, ( a ) and (b), 
and the angle (/f) as defined in Pig. 10. 

For the angle 0, it can be shown in a similar manner as above that: 

R - (ab)/ \/b 2 cos 9 - a 2 sin 8 ( 29 ) 

The development of the equations for the first and second derivatives 
of Eqs. (28) and (29) shall not be attempted here, as they become too complex 
• to be practical for layout work. When the need arises for such expressions, it 
is simpler to convert to the rectangular coordinate system, or to the use of 
power curves which will be outlined in the next chapter. 

1 - 40 - Other Basic rTfiomfttTH r- shapop.- Other basic geometric 'shapes, such as 
the parabola, hyperbola, etc., could also be used -in developing contour 
lines for various shapes of bodies. However, by using the power curves of the 
next chapter, the circle and the ellipse can be so easily modified to fit any 
desired contour that there is no need, at this time, to consider any other 
basic geometric figures. Later on a discussion is given of some applications 
of the parabola, which may be used sometimes to advantage in the development 
of equations. 
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CHAPTER II 
THE POWER CORVES 

2.10 - 'Tjyhroductiort s- Power curves are so called because the powers, or 
exponents, of the dependent variables of the equations are the main factors 
that govern the shape of the curves. The equations may be single or multi-term, 
and each consists essentially of four parts; the dependent variable, the 
independent variable, the constant, and the exponent. TOien either of the latter 
two is unity it can be, and usually is, omitted in the equation. Any symbols 
may be used for the component parts, but to follow the general practice of 
mathematical study, it shall be assumed that (y) is the dependent variable, 
(x) the independent variable, (k) the constant, and (n) the exponent, or power. 
Thus a single term power equation will be written: 

y-k(x) n .............(30) 

The dependent variable (y) is a function of (x), and can be computed 
for any (x) desired. It represents the ordinate, while (x) is the abscissa. 

The constant (k) governs the height, or thickness, of the curve, and 
is usually small when the exponent is large, and large when the power is small. 
A change in magnitude of the constant will increase or decrease all (y f s) in 
the same proportion. « 

The exponent, (n), which governs the curvature, may be divided into 
three classes. First, the exponent can be equal to unity. A term with this 
power will form a straight line, the constant of which determines the slope. 
Second, the exponent can be larger than unity. A term with this power will form 
a curve that follows the X axis at the origin for an infinitesimal distance, 
turns up swiftly, and continues with ever-decreasing curvature. Third, the 
exponent may be smaller than unity. A term with this power will form a curve 



18. 

that follows the Y axis at the origin for an infinitesimal distance, turns 
swiftly to the right, and then also continues with ever-decreasing curvature. 
This exponent is usually required for the first term of all equations of 
streamlined shapes developed by this method. 

If Eq. (30) is evaluated, it will be found that, except at the origin, 
the slope of the curve will never be zero or unity, no matter how far it is 
extended. The curve, therefore, will never be parallel to either the X or T 
axes. To obtain this, however, one or more terms may be added or subtracted. 
Thus, if a negative term with a small constant and an exponent larger than 
unity, (if that of the first term is less than unity), is added, the curve will 
become parallel to the X axis at some point and will then eventually intersect 
it. If a positive term is added under the same conditions, however, the curve 
will not become parallel to the T axis, as will be shown in a later chapter. 

An equation with one or more termg added to the first is called a 
multi-term equation. If subscripts corresponding to the term number are used 
for (k) and (n), a multi-term equation will then be: 

y - fc^x}* 1 - k 2 (x) n 2 - k 3 (x) n 3 ...etc (31) 

The slope at any point on the curve can readily be determined by 
means of the first derivative, which is: 

(dy/dx) - n^x)" 1 ^ - n^tx)" 2 " - n 5 k 5 (x) 5 ...etc ...(32) 

It can be easily seen that Eq. (32) is in the same form as Eq. (51), 
since the values (n-^), (ngkg), etc., are also constants. Similarly, the 
equation for the second derivative is expressed in the same form. These equation 
forms, while they may seem involved at first glance, are simple and have an 
unparalleled flexibility ideal for layout work. 
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It has been found that it is easier to deal with power curve equations 
with fractional exponents when common logarithms are employed. Therefore, in 
the discussions to follow on the use of power curves, the constants are left 
in log form and noted with the letter (N), signifying that the number follow- 
ing is not the actual value of the constant but the logarithm of the constant. 
For example, a single term equation with a constant of 5 and an exponent of 

0.5 would be written: 

0.5 
y = N(0.47712)(x) 

For evaluation of (y) in the above equation, the log of (x) is multi- 
plied by the exponent 0.5 and added to 0.47712. The antilog of the sum, as 
found from logarithmic tables, then equals the value of (y). 

In the case of an exponent of unity, there is no need to resort to the 
use of logarithms, and the actual values of the constants can be more readily 

used. - ' 

2 .20 - »hg cm™i« *nd » j.o gmn»» it, Power ffrr™ Form:- The use of power curves, 
and a better understanding of their form, can best be shown by first working 
with the equations of the circle and the ellipse as derived in Chapter I. 

Consider the equation of the ellipse with the origin on the X ax?, s and 
on the periphery, as given by Eq. (25). This can be written: 

. r / » (b 2 /a*)(2ax-x 2 ) 
or, multiplying out the right hand side: 

/ = (8b 8 x)/(a) - CbV)/(«*> 
The right hand side of this equation can now be put in the form of the 
right hand side of Eq. (51), or in power curve form, by letting: 

■ ^ . (»*)/(a) • - (5M 

lfc * (bW) ( 52b > 
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Then the equation may be written: 



y 2 = k^x) - k^OO 2 
or, solving for (7): 

7 » (^x - k^x 2 ) 2 (35) 



Thtis Eq. (35) represents the power curve form of Eq. (25). In this I 



form, as will be seen later, the equation becomes more flexible for use in 
modifying the basic elliptical shape into contour lines of different curvatures. 

The slope at any point on the curve of this equation can be found, as 
before, by means of the first derivative. Hence: 

(dy)/(dx) - Kkjxr - kgx 2 )" 2 ^ - 21^0 
_ (ki - 2kpx) 
2(k,x - kj>x 2 ) 2 
Substituting (y) for (k^x - k^c 2 ) 2 , from Eq. (33): 

(dy)/(dx) - (k^ -2^/(27) (34) 

By substituting in Eq. (34) the values of (kj_) and (kg), as given 
above, it will be found that it is the same equation as Eq. (26), but that it 
is in a somewhat simpler form. 

The second derivative is found as follows: 

(dy)/(dx) = Kkjx - kgxV^ - 2kgx) 
Then^ differentiating the above: 

d 2 y/dx 2 = i [(fcjx - k 2 x 2 )-2(-2k 2 )] + ^([k! - 2kgx][^Cklx - k^ 2 )" 5 /^^ - 2k«>x)] J 



-2kg (k x - 2kgx) 



2 



2(k 1 x - kgx 2 ) 2 4 (k 1 x - kgx 2 ) 3/2 

_ 4kg (k x x - kgx 2 ) + (k-L - 2k g x) 2 
4(k x x - kgx 2 ) 5/2 



>i 
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- 4k g (k 1 x - k 2 x 2 ) -i- (Ig - 2k 2 x) 
4(k 1 -k 2 x 2 )(k 1 x-k 2 x 2 )2 

Substituting (y) for (k x x - k«>x 2 )^, from Eq. (S3), and (/) for 
(k.x - k p x 2 ), the final equation in power curve form becomes: 

-4k 2 y 2 + (k ^ - 2kgx) 2 

(dV^ 2 ) » 73" 

4y 



.(55) 



2.21 - Hombining Part of an Ellipse wi th an Arc of a Circle:- As an 
example of the way in which the power curve form of the equation of an ellipse 
is used, consider the case where it is necessary to add part of an ellipse to 
the arc of a circle. This can be done mathematically for layout work, and should 
be done in this manner to assure a smooth continous curve, with a uniformly 
changing rate of curvature * 

In Fig. 11, assume that the radius of the circle is (r), the distance 

D. 



t = slope 




FIGURE 11 
from the center line to the maximum width required is (c), the distance from 
the periphery of the circle to the point of tangency is (d), and the distance 
from the maximum width desired to the point of tangency is (b) . 
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Eq. (7aX gives the equation of a circle with the origin at (B), as 
in Jig. 11. 

y=(2rx-x 2 ) 2 (7a) 

The slope at any point is given by Eq. (9): 

(dy/dx) » (r - x)/y ........(9) 

Using Eq. (7a), the value of (y) at the point of tangency, (D), can be 
found, and Eq. (9) permits the slope to be found at the same point. Having 
found these values, the equation for an ellipse passing through point (D) with 
the given values of (y) and (dy/dx) can now be developed. Using preferably the 
power curve form of the equations for an ellipse with the origin at (A), as in 
Fig. 11, the right hand side of Eq. (35) is equated to the right hand side of 
Eq. (7a), and also the right hand sides of Eqs. (54) and (9) can be equated. 
Hence,, for the ordinate at the point of tangency: 

(2rx - x 2 )^ = (^x - kj>x 2 ) 2 
or, using the symbols as in Fig. 11, and cancelling the root signs: 

(grd - d 2 ) = (kjb - kgb 2 ) . (36) 

Similarly for the slope: 

(r-x)/y * (1^ -2kgx)/2y 
or: 

2(r - d) = (k x - 2kj,b) (37) 

Eqs. (36) and (57) can now be solved simultaneously to determine the 
■values of k x and kg, and these values then put in the form of Eq. (33) to give 
the required equation of the ellipse. 

An example is given below to illustrate the method outlined. Assume, in 
Fig. 11, that the following quantities are known or given: 

b = 35 
c = 40 
d = 25 
r - 50 

Using Eq. (36): 

2(30)^25) - (25) 2 - ^(55) - kg^) 2 



,im 
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or: 

25 - k x -35(kg) (38) 

Using Eq. (37): 

, 2(30 - 25) - k^ - 2kg(35) 
or: 

'-■ K> -kjL - TOOcj,) (39) 

Solving Eqs. (38) and (39) simultaneously for k-, and k« by first 

multiplying Eq. (39) by 2, we get: 

50 -2(1^) - 70(kg) 

10 - (kj) - 70(kg) 

40 - (k x ) 

Substituting this value of k. in Eq. (58), (or Eq. (59) could also be 
used), the value of kg is found to be: . 

25 m 40 - 35(kg) 

kg = 15/35 
= 0.42857 
Hence, the aquation of the required ellipse, in power curve form, is: 

y = (40b - 0.42857^)2 
or, since (b) really represents the conventional (x) term: 

y = (40x - 0.42857x* : ) 2 (40) 

This is then the equation of an ellipse with the origin on the X axis 
and on the periphery, (point (A) in Fig. 11), which is also tangent with the 
periphery of the given circle at point (D). The slope at any point on this 
ellipse is then obviously the first derivative of Eq. (40). 

2.22 - Combining Two Ellipses :- Many times it may be necessary to 
enlarge the cross section of elliptical bodies, with the least amount of rework 
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of existing parts and structure. Such may be the case with an airplane fuselage, 
for example, when more space is required without changing the major fuselage 
longerons, or affecting too much of the internal layout. 

The simplest way of doing this is to add a semi-ellipse, with a larger 
major axis, to the ellipse on hand. In Pig. 12, below, let it be assumed that 
it is desired to add a portion of an ellipse to the original elliptical cross- 
section so as to increase the height by 10 inches. (For the sake of maintaining 
the same equations as previously developed, the "height" in this case is measur- 
ed horizontally along the X axis. Obviously, for such layout work, it is 
immaterial in what direction the "height" is considered). 




FIGURE 12 
The known, or given, dimensions of the original and desired cross- 
sections are indicated on Fig. 12. Since the (y) values and slopes of both 
ellipses must be the same at Station (B) In order for them to be tangent, it 
is only necessary to equate the respective "curve" and "slope" equations for 
the two ellipses. Therefore, letting ki and lc^ be the constants for the original 
ellipse, and k, and kp for the segment to be added, we have, from Eq. (53) s 
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(k£x - kg* 2 )* = (kj'x - kg**¥ (41) 

and, from Eq. (34): 

(kj_ - 2kgx)/2y ~ d4 f " 2*4'x)/27 "' < 42) 

It must be remembered that the values of (x) are not the same when 
used for the original and new ellipses, since the original ellipse has its 
origin at (C), and the new segment to be added has its origin at (D). 

The constants of the original ellipse can be found from Eqs. (52a) and 

(52b): 

k{ - (2b 2 )/a = 2(50) 2 /40 = 45 

k! - (b 2 )/(a 2 ) - (30)?(40) 2 = 0.5625 
2 

Substituting these in Eq. (41), and remembering that the value of (x) 
to the point of tangency for the original ellipse is 28 inches: 
[45(28) - 0.5625(28) 2 ]2 = (38k^' - 1444kg') 2 



or: 



k[' - 58kg' = 21.5526 (43) 



Substituting the values of k-j^ and kg in Eq. (42): 
45 - 2(0.5625)1,28) - k^' - 2(38)kg 



or: 

it it 



1^ - 76kg = 13.5 (44) 



The value of kg 1 can now be found by solving Eqs. (45) and (44) simul- 



taneously, giving: 



kg* s 0.2119 



T t 



k{ can now be found by substituting this value of kg in either Eq. 
(45) or (44): 



k^'= 29.604 
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The power curve equation of the new segment of the ellipse, with its 
origin at (D) and tangent to the original ellipse at Station (B), is then: 

7 - (29.604* - 0.2119x E )2 

2.30 - TtovftlnnraftT it. nf a Mr.crTp Tm-m Pn« 01 > fhmrp TCqnp-H ™«- The previous 

sections have shown how the power curve equations can be used in the case 
where it is required to pass the curve of an ellipse through a point, the 
ordinate and slope of which are given. It shall now he shown how the equation 
for any type of a curve can be developed which will pass through one or more 
given points, and, at first, the single term equation will be considered. 

The standard single term power curve equation has already been shown in 
Eq. (30) to be of the form: 

y = k(x) n ... ..(30) 

The parabola is an excellent example of the above form of equation. 

The standard parabolic equation is: 

y 2 = 2px (45) 

where (p) equals twice the distance from the focus to the vertex, 
or origin. The above equation may be written: 

y = "(Ep)*(x)* 

A 
and, if (k) be substituted for (2p) 2 ; 

y = k(x)°' 5 (46) 

This is a single term power curve equation of the same form as Eq. 
(30), and is often used, as is, for the first term of the multi-term equations 
to be developed later. 

Referring to Fig. 13, let it be assumed that it is required to pass 
a curve through the two points (P) and (Q), whose ordinates and abscissae are 
as shown. 
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y FIGURE 15 
If it first be assumed that a separate equation, in the form of Eq. 
(30), be written for point (P) and also for point (Q), these equations become: 

b' akUV ......(47) 



B'^kU")* 



.(48) 



Then, from Eq. (47): 



k- (B')/U f ) n - < 49 > 



and, from Eq. (48): 



k S (B M )/(A ,, ) n 



.(50) 



Since it is now desired to combine Eqs. (47) and (48) into a single 
equation, there can be only one constant, or, in other words, the values of 
(k) in both equations must be one and the same. We may then write: 



or: 



.It /_T, 



(A ,, /A , ) n = (B'VB 1 ) 
Transposing, and using logarithms, the value of (n) is found to be: 

n - logtB'VB'Vlog^'A') .....(51) 

If this value of (n) is now substituted back in either Eq. (49) or 
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(50), the value of (k) can be found, and then the equation is developed in 
the familiar form: 

y = k(x) n 

Often, however, it happens that a second point on the required curve 
is not known, or is of no importance, but instead the slope of the curve at a 
point is known. Then the derivative of the equation is employed for this 
point. For example, in Fig. 14, let (A) be the abscissa and (B) the ordinate 
of the point (P), and let (t) be the slope at this point. 

T P^—-*'^ t = slope 




y FIGURE 14 
Substituting in the standard Eq. (50) these values, we get: 



B = k(A) n 



or: 



The first derivative of Eq. (52) is obviously: 

(dy/dx) = (dB/dA) » t = nk(A) n " 1 
The value of (k) is given by Eq. (53), and hence: 
t = n(B/A n )(A) n _1 

= nB/A 
and then: 

n = At/B 



(52) 



k = B/(A) n (53) 






.(54) 
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Thus, knowing the values of (k) and (n) fro., fts. (55) and (=4), 
respectively, the standard alngla tan. power carve aquation can be set up. 

Should the Slope, (t), be of no importance, then the exponent (n) can 
be assumed and (k) tnd from ft. (55) a*d both values used in ft. (50) to 

develop the curve. 

fts. (49) or (50), (51), (55), and (54) for the Boat part form the 
b asis for the development of power curve equations, and should be memorised. 

, 51 - .,-iW.i- of a «rr'--~ —- "■"* *"**"" " *" ^ 

* *« . the u»« of th. equations developed in the preceeding 
purpose of demonstrating the use ol tne equa., 

section, let it be assumed that it is required to form an equation for a curve 
passing through two given points. Therefore, in Fig. 15, assume the following 
to be given: 



A - 



15 
10 



A ,fl =32 
B " =14 
nrst substitute in Eq. (51) to find the exponent needed: 
n - log(14/l0)/log(52/15) 

I n F ( 14> - log (10) 

" log(32) - log (15) 

= (1.14615 - 1.000)/(1.50515 -1.17609) 

. 0.44408 
The value of 00 is «— from ft. (55); the single primed values will 
b e used, although the double primed ones could be used just as well: 

i-« , \//. , \0.44408 
k = (B )/(A ) 

or, using logarithms for ease of computation: 

k-« log(B') - 0.44408(log k) 

= log (10) - 0.44408 (log 15) 

s 1.0000 - 0.44408(1.1761) 

= 0.47772 
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The required equation i3 then: 

y=N(0.47772)(x) ' 44408 
As has been noted previously, the reader should be careful to remember 
that, in the above equation, the letter (N) signifies that the real value of 
the constant is not 0.47772, but is the number whose logarithm is 0.47772. The 
only reason the equation is left in this form is to facilitate computation of 
various values of (y) for assumed values of (x). For example, if (x) « 10, then: 

y - N(0.47772) + 0.44408(log 10) 
- N(0.47772) 4 0.44408(1) 
■ N(0.9218) 

The number whose logarithm is 0.9218 is found from log tables to be II 

8.35, and hence: t| 

y s 8.35 when x = 10 ft 

'"•;;'■■ ■ . ■• ■ ■ ' ' ".' : | 

As a second example, consider the case where it is desired to find the | 

- " : - ^ 

equation of a curve which passes through a point at a given slope, (t). Referr- 
ing to Fig. 14, let: J 



From Eq. (54): 



A = 20 
B = 10 
t = 0.2 



n ■ At/B 

= (20)(0.2)/lO 

■ 0.4 



and, from Eq. (53): 



k - B/(A) n 



= N §-og(lO) - 0.41og(20)] 
= N(0.47959) 
Hence, the required equation is: 



rf : ; 



=10/(20) - 4 I 
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B = N(0.47959)(A)°' 40 
or, in usual terms: 

The equation for the curve can easily be checked for any error as 
follows. Substituting A * 20, then: 

B » N(0.47959)(20)°* 4 ° 
:= N(0.47959 + 0.40 log (20)) 
-- '-" « N(0.47959 + 0.52041) 
» N(1.0000) 
-'- io. 000, which checks. 
To check the slope, consider the first derivative of the equation: 

(dy/dx) j* (dB/dA) - t = (0.4) [ff(0.47959)](x)~°* 

■ ' - N(log(0.4) +■ 0.47959 - 0.6 log (20)) 
= N [(0.60206-1) + 0.47959 - 0.6(1.30103)] 
- N(0.30103-l) 
= 0.2, which checks. 
With variable slopes, the number of curves possible through two points, 
or, as in this case, the origin and a point, is infinite. For each assumed 
exponent, (n), the constant can be solved for by means of Eq. (53). To show, 
and get acquainted with, the general trend of these curves, a series of equations 
has been developed and the ordinates computed and plotted. 

In Table II, the equations are given for curves passing through the 
origin and the point A - 40, B - 40, for different assumed values of the 
exponent, (n). The values of (k) were found from Eq. (53), and then put in log 

form. 
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TABLE II 
Curves Through Origin and Point (A = 40, B « 40) 



Assumed 
Exponent 


k 


Equation 


0.3 


N(1.12144) 


y = N(1.12144)(x)°* 3 


0.5 


N(0.80103) 


y = N(0.80103)(x)°* 5 


0.7 


N(0.48062) 


y = N(0. 48062) (x)°* 7 


0.9 


N(0. 16021) 


y = N(0.16021)(x)°' 9 


1.0 


N(o. 00000) 


y = x 


2.0 


N(0. 39794-2) 


y ="n(0 ; 39?°4-2)(x) 2 


4 -° 


N(0.19382-5) 


y - N(0.19582-5)(x) 4 



Table III shows the solution for values of (y), in the equations 
given in Table H, for assumed values of (x). These equations are then plotted 
on Fig. 15. 

TABLE III 
Curves Through Origin and Point (A = 40, B = 40) 



'I 



-X.. 


(y) for an exponent (n) of: 


„Pt3., 


0.5 


0.7 


0.9 


1.0 


2.0- 


4.0 


1 


13.227 


6.325 


5.024 


1.446 


1.000 


0.025 


0.000 


2 


16.284 


8.944 


3.024 


2.699 


2.000 


0.100 


0.000 


5 


21.436 


14.142 


9.330 


6.156 


,3.000 


0.625 


0.010 


10 


26.390 


20.000 


15.157 


11.487 


10.000 


2.500 


0.156 


20 


32.490 


28.285 


24.623 


21.436 


20.000 


10.000 


2.500 


30 


36.693 


34.641 


32.704 


30.876 


30.000 


22.500 


12.656 


40 


40.000 


40.000 


40.000 


40.000 


40.000 


40.000 


40.000 



4 I 
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FIGDEE 15 



For the second example, equations are developed for curves drawn 
through the origin and the point A = 100, B = 40. Table IV shows the equations 
developed in the same manner as in the previous example. 

The ordinates for the equations in Table IV are given in Table V, and 
the results are plotted on Fig. 16. A comparison with Fig. 15 shows clearly 
that, for the same exponent, the curves have the same general trend, those of 
Fig. 16 being flatter due to the greater (A) distance. The reader will find it 
time well spent to study these curves and become familiar with the way in which 
the shape of the curves of a single term power curve equation may be modified 
by using various values of the exponent. 
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TABLE IV 
Curves Through Origin and Point (A = 100, B = 40) 



Assumed 
Exponent 


k 


Equation 


0.5 


N(1.00g06) 


7 = N(l.00206)(x)°* 3 


0.5 


N(0. 60206) 


y = N(0.60206)(x) * 


0.7 


N(0.20206) 


7 * N(0.20206)(x)°* ? 


0.9 


J N(0.80206-l) 


7 = N(0.80206-l)(x)°* 9 


1.0 


N(0.60206-l) 


7 = N(0.60206-l)(x) 


2.0 


N(0.60206-3) 


7 - H(0. 60206-3) (x) 2 


4.0 


N(0. 60206-7) 


7 = N(0.60206-7) (x) 4 



TABLE V 
Curves Through Origin and Point (A = 100, B = 40) 



X 


(7) for an exponent (n) of: 


0.3 


0.5 


0.7 


0.9 


1.0 


2.0 


4.0 


5 


16.280 


8.940 


4.910 


2.700 


2.000 


0.100 


0.000 


10- 


20.050 


12.650 


7.980 


5.040 


4.000 


0.400 


0.004 


• 20 


24.680 


17.890 


12.970 


9.400 


8.000 


1.600 


0.064 


40 


30.390 


25.300 


21.060 


17.530 


16.000 


6.400 


1.024 


60 


34.320 


30.980 


27.980 


25.260 


24.000 


14.400 


5.184 


80 


37.410 


35.780 


54.220 


32.720 


32.000 


25.600 


16.390 


100 


40.000 


40.000 


40.000 


40.000 


40.000 


40.000 


40.000 



1 
I 



H 
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2.40 - T wo Term Power Curve Equations : - It has been stated previously that a 
Bingle term power curve can not be made to parallel either the X or Y axis, 
(except for a minute distance at the origin), and therefore can not be made to 
intersect either of these axes. 

To obtain a curve that will parallel either the X or Y axis at some 
point, and later intersect the axis, one or more terms must be added to the 
first term. The equation is then called a multi-term equation. The constants 
of all these terms can be solved simultaneously, while the required exponents 
either have to be assumed or found by a trial and error method. 

The two term equations are by far the simplest of the multi-term 
type, as they require little work not only in their development but also in 
the computation of the ordinates. However, two term equations can only be 
directly developed when certain information is given, and, when this informat- 
ion is not available, resort must be made to the trial and error method out- 
lined in Section 2.60. 

An analysis of two term equations shows that the exponents are the 
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only factors that govern the location of the maximum ordinate, and that the 
constants control the thickness ratio, where thickness ratio is defined as the 
ratio of the maximum ordinate to the distance from the origin to the point of 
intersection with the axis. 

The following sections demonstrate the manner in which two term equat- 
ions can be developed directly. 

o AT - n»+.»wm-Ti«t1on o f t.hft Const p nt.s of a Two T*rm PoWfiT CUTTS 
fr p i pMnn When t.he Kxponer rhs are Known; - The simplest condition for which a 
two term power curve equation can be written is when both exponents are given, 
or assumed. Such conditions arise in practice through experience with the 
type of curves obtained with various exponents, it very often being desired to 
duplicate the general shape of a curve previously developed but instead pass- 
ing it through two different points. 

However, only two points on the desired curve can be handled in develop- 
ing the two term equation, since a requirement for passing the curve through 
three or more points necessitates the use of three or more term equations. 

The constants of an equation of two terms can be solved simultaneously 
in the same manner as the exponent and constant of a single term equation. It 
will be demonstrated here, by means of a numerical example, a method of solving 
for the constants of a two term equation, the exponents being assumed as 0.5 
and 1.0 for the first and second terms, respectively. The curve is required to 
start at the origin and pass through the points (x) - 25, (y) - 12, and (x) = 
100, (y) - 0. This can be assumed to be a typical example of developing the 
equation for an entire streamline shape which starts at the X axis and inter- 
sects the X axis again at a point 100 units away. 

Eq. (31) shows the standard form of the multi-term power curve 
equation, and therefore, substituting the coordinates of the two points. to 
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form the two preliminary equations, we get: 

y = k^xfl - k 2 (x) D 2 



.(31) 



and: 



12 = ^(SS) * 5 - kgCSS) 1 * 



o = k^ioo) * 5 - kgCioo) 1 * 



Solving these two equations simultaneously by first extracting the 
roots and multiplying the first equation by 2 and subtracting it from the 
first, (kg) is found to be: 

k 2 = 0.48 

= N(0.68124-l) 
The value of (k x ) is found by substituting this value of (kg) back in 

either one of the above equations: 

K - 4.80 

= N(0.68124) 
The desired equation is then, from Eq. (31): 

y = N(0.68124)(x)°* 5 - N(0.68124-l)(x) 

2.42 - Ttet.P-nrrmati on of the fr monents of a Two Term Power Curve 
pquation : - Two term power curve equations can be readily written for the 
condition where the value and location of the maximum ordinate , and the length 
along the X axis from the origin to the point of intersection, are given. This 
condition often occurs, as it actually consists of developing a streamline 
shape of a given length, thickness ratio, and position of maximum thickness. 

In this case, curves are developed in this section to permit the ex- 
ponent of the second term to be found directly, after the first term exponent 
has been assumed from a consideration of the type of curve desired. (See Figs. 
15 and 16). 
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As mentioned in Section 2.40, the exponents of a multi-term power 
curve equation govern the location of the maximum ordinate. The location of 
the maximum ordinate along the axis is usually expressed as a percent of the 
distance from the origin to the point of intersection of the curve with the 

axis. Thus, in airplane airfoil section layout work, for example, the distance '% 

from the origin to the intersection of the curve with the X axis is called the J 

"wing chord", and the maximum ordinate, or airfoil thickness, is located at | 

-■■-■■■■-;-■■ :■■■■-■■■ } 

some percentage of this chord from the origin, or leading edge. 1 

To show how the exponents of the terms control the location of the | 

1 

maximum ordinate, consider the following two term power curve equation, with | 

the constants equal to unity: 1 

y=(x)°- 5 -(x) (55) 1 

From calculus, it is known that if the first derivative of an equation S 

is set equal to zero, and solved, the value of the dependent variable, (y in | 

this case), is a maximum or minimum. Hence, taking the derivative of Eq. (55): J 

(dy/dx) = 0.5(x)"°* 5 - 1 I 

and, letting this equal zero, the value of (x) becomes: 

- .-0.5 
0.5(x) -1=0 

'"■'... (0.5)/(x)°* 5 = 1 } 

x = 0.25 

This is the location along the X axis at which the ordinate, (y), is .( 

■ •"■"• . ■■ ■ -'I 

a maximum. The point on the X axis where the curve intersects the axis is J: 

h -." _'■?#:£. 

obviously where (y) is equal to zero, and therefore, setting Eq. (55) equal to i 

• 'I 

zero and solving for the value of (x) at the intersection: • » 

■ ■ ■ ■ i 

(x) ' 5 - x = J 



•# 
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(x) - occurs at the origin, while the other point at which (y) is 
zero is at a point one unit away from the origin. The maximum ordinate is 0.25 
units away from the origin, or therefore 25^ of the distance from the origin 
to the point of intersection. 

The constants in Eq. (55) were assumed to be equal to unity, but a 
change in constants will not affect the location of the maximum ordinate, 
assuming the exponents to be kept the same. This is shown in Table VI, where 
the same exponents as in Eq. (55) are used, but the constants changed. In 
this Table, the term "chord" is used to designate the distance from the origin 
to the point of intersection of the curve with the X axis. 

TABLE VI 



Effect of Constants on Location 


of Maximum 


Ordinate 


Equation 


Location of 

Max. Ordinate 

(Ifaits) 


Chord 
(Units) 


Location of 
Max. Ordinate 

(% Chord) 


y - (x)°- 5 - 2(x) 


0.0625 


0.25 


25 ^ 


y = 2(x) - 5 - 2(x) 


0.2500 


1.00 


25 


y = (x) * 5 - 0.5(x) ; 


1.0000 


4.00 


25 


y = (x)°* 5 - 0.2(x) 


6,2500 

■ 


25.00 

• 


25 



The computations necessary for completing Table VI were made in 
exactly the same manner as outlined for Eq. (55). 

Table VII shows, however, that when the exponents are changed, the 
location of the maximum ordinate is also changed. It will be noted in Table 
VII that the exponents of the fourth and eighth, and the sixth and ninth, 
equations are the same, and therefore the locations of the maximum ordinate 
in percent of the chord are the same. 
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TABLE VII 
Effect of Exponents on Location of Maximum Ordinate 



Equation 


Location of 

Max. Ordinate 

(Units) 


Chord 
(Units) 


Location of 

Max. Ordinate 

(% Chord) 


, ,0.5 , .1.5 
y - (x) - (x) 


0.3333 


1.000 


33.33 


, x0.5 . x2 
7 = (x) - (x) 


0.3968 


1.000 


39.68 


7 - M°' 5 - (x) 2 * 5 


0.4484 


1.000 


44.84 


y = (x)°- 5 - (x) 3 


0.4883 


1.000 


48.85 


/ \0.*3 / \3. o 
7 - (x) - (x) 


0.5227 


1.000 


52.27 


/ A 0.5 / v4 
7 - (x) - (x) 


0.5521 


1.000 


55.21 


, .0.5 . .10 
7 - (x) - (x) 


0.7296 


1.000 


72.96 


, x 0.5 . .3 
y - 2(x) - 0.5(x) 


0.8505 


1.741 


48.83 


7 - 3(x)°* 5 - 0.0000006(x) 4 


45.2900 


82.023 


55.21 



It should be remembered that if the constants of two or more equations 
are equal, the thickness ratios, (as defined in Section 2.40), are the same, 
regardless of whether or not the exponents are the same. Hence, in the first 
seven equations of Table VII, the thickness ratios are equal. 

In the equations in Table VII, only the exponents of the second terms 
were changed. However, it Is also true that the location of the maximum 
ordinate will vary if the exponents of the first terms are changed while the 
exponents of the second terms are unchanged. A general two term power curve 
equation will now be developed for this case where the value and location of 
the maximum ordinate , and the distance to the axis intersection, are given. 
The distance from the origin to the point where the curve intersects the axis 
will be called the "chord", and, if for the sake of consistency with the rest 

of the text the axis is assumed to be the X axis, then (y) = at the 
intersection. 
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The maximum ordinate, or the maximum value of (y), occurs where 
the curve is parallel to the X axis, or, in other words, where the slope, 
(dy/dx), is equal to zero. Thus, utilizing these definitions, two equations 
can be written giving the relationship of the exponents to determine the 
location of the maximum ordinate along the chord. Using the general form of 
a two term equation from Eq. (31), then, and writing it for the point at the 
end of the chord where (y) is equal to zero, we have for the first equation: 

y = k^xA - k 2 (xf2 = ......."(56) 

The first derivative of Eq. (56) gives the slope, which is equal to 
zero when (y) is a maximum. Hence, for the second equation: 

(dy/dx) = n^x)* 1 " - n^x)* 2 =0 (57) 

If (L), for length of the chord, is substituted for (x) in the first 
equation, and (CL) is substituted for (x) in the second equation, where (C) 
represents a constant whose value is the position of the maximum ordinate in 
hundredths of the chord, we get for Eq. (56): 

^(L) 11 ! = k 2 (L) n 2 
or: 

(fcjAj,) = (L) n2 /(L) ni (58) 

Eq. (57) becomes: 

- 1 n„ - 1 



n-jk-^CL) 111 " = n 2 k 2 (CL) 2 



or: 
(k^kg) = n^CL)* 2 " 1 /n 1 (CL) il " (59) 

The relationship between the exponents can now be found by equating 
Eqs. (58) and (59), as follows: 
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(L)"* n 2 (CL) n 2- 



n«-l 



(L) n l n n (CL) n l _1 



(L) n l (L) D 2 



b-lCC) 11 !- 1 = n^C)^- 1 



or: 



(nj/njj)-. = (C) D 2 n i ....(60) 

Eq. (60) readily permits the determination of the location of the 
maximum ordinate for any two term power curve equation if the exponents are 
known. For example, if (bt) and (n«) are assumed to be 0.5 and 1.0, respect- 
ively, the location of the maximum ordinate is found as follows: 

(o.sA.0) = (c) 1 - - - 5 

0.5 = (C)°- 5 

log (0.5) = 0.5 log (0) 

log (C) = (0.69897-l)/(0.5) 

= 1.39794-2 

or : 

C - N(0. 39794-1) 

» 0.25 

Hence the maximum ordinate is located at (x) = 0.25(c), or at 25^ of 
the chord from the origin. 

Further evaluating Eq. (60) for fixed values of (n x ) and various values 
of (n«), in the same manner as in the above illustrative example, a series of 
curves can be computed and plotted from which (n 2 ) can be easily found for any 
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maximum ordinate location desired. This has been done in Fig. 17. 

For demonstrating the method of developing two term power curve equat- 
ions by the use of Fig. 17, a few numerical examples will be given. First, let 
it be assumed that the equation to be developed is to have a first term 
exponent of 0.4, and that the curve is to have a chord of 262 inches, a maximum 
ordinate of 26 inches, and the location of the maximum ordinate is to be at 
(x) = 112 inches. 

Then the may*™™ ordinate is at: 

C - (112/262) (L) 
= 0.4275 (L) 

or: 

= 42.7556 of the chord 
From Fig. 17, for (n x ) = 0.4 and 100(C) = 42.75, the required value of 

(n 2 ) is: 

n 2 = 2.625 

The equation is now: 

y = k-^x) - k 2 (x) 

All that now remains is to solve for (k x ) and (kg), which can be done 
ob the manner outlined in Section 2.41. Hence, at the location of the maximum 
ordinate, y = 26 inches and x - 112 inches, Eq. (61) becomes: 

y = k l( 112) - 4 - k,^) 2 - 625 . 26 • («) 

The second equation is written for the point of intersection of the 
curve with the X axis, where (y) - inches and (x) = 262 inches. Therefore: 

y = ^(282)°-* - k^) 2 ' 625 = («> 

Eqs. (62) and (65) become, after raising to the required powers and 
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extracting the roots, (which, as mentioned previously, can usually be done 
most easily by the use of logarithms); 

y - 6.602^) - 239442(kj>) - 26 

and: 

y - 9.2752(1^) - 2228632(kg) = 
Solving these equations simultaneously, (k x ) and (kg) are found to be: 

1^ = 4.63834 
= N(0.66656) 

and: 

k = 0.000019304 

= N(0.28565-5) 

Eq. (61) can now be put in terms of the final equation as follows: 

n A „»# x2.625 f 64 \ 

y = N(0.66636)(x)°' 4 - N(0.28565-5)(x) I**' 

In a Similar manner as shown above, the use of a value of (b^ = 0.5 
instead of 0.4, with the chord, maximum ordinate, and location of the maximum 
ordinate the same as in the first example, produces the following equation: 

y -. H(0.49636)(x)°- 5 - N(0.14341-4) (x) 2 ' 3 — < 65 > 

If (^ is taken as 0.6, and all other variables are again assumed the 
same as in the first example, the equation becomes: 

y m N(0.33520)(x)°* 6 - N(0.82865-4)(x) 2 * -( 66 > 

Table Till presents an evaluation of Eqs. (64), (65), and (66) for 
various assumed values of (x), and the values from this Table are plotted on 
Fig. 18. 



Ub. 



TABLE VIII 
Evaluation of Eqs. (64), (65), and (66) 



(*) 




(7) 






Eq.(64) 


Eg. (65) 


Eq.(66) 


1 


4.64 


5.14 


2.11 


2 


6.12 


4.45 


5.50 


5 


8.83 


7.01 


5.67 


10 


11.64 


9.89 


8.54 


25 


16.72 


15.45 


14.43 


60 


22.96 


22.58 


22.26 


112 


26.00 


26.00 


26.00 


162 


23.52 


25.11 


22.99 


204 


16.61 


16.24 


16.01 . 


238 


7.97 


7.68 


7.51 


262 


0.00 


0.00 


0.00 




7 - N(0.66636)(x)°* 4 - N(0.28565-5)(x) 2 '^l 
7 = N(0.49635)(x) 0,5 - N(0.14341-4)(x) 2 * 3 2- 
= N(0. 33520) (x) * 6 - N(0. 82865-4) (x) 2 *°£ 




10 25 



112 



162 



204 



238 262 



FIGURE 18 
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By comparing the curves of the three equations, it becomes at once 
apparent that an equation that has a large exponent for the first term gives 
a curve with narrow, or more pointed, front and aft sections, and that the 
bluntness increases as the exponent of the first term decreases. In all these 
equations the first term could be placed after the second, thus making the 
second term the first, but the equation would still be the same. In this work, 
however, it shall be understood that by the first term is meant that term 
which has the lowest exponent of (x), which, in most cases, is less than unity, 
and that the term number is designated in the order of the increase of the 
exponents of (x) . 

The illustrative examples presented above show the means whereby two 
term power curve equations can be developed when the maximum ordinate conditions, 
and the chord, are specified. It has also been demonstrated how the shape of 
the curve can be modified by the choice of the exponents, and by the values 
chosen for the maximum ordinate and its location. 

2.43 - Other Conditions for which Two Term Equations Can be 
Developed ;— It has been shown in Sections 2.41 and 2.42, respectively, the 
manner in which two term power curve equations can be written when (1) the 
exponents are given, or assumed, and (2) the magnitude and location of the 
maximum ordinate, and the chord, are given. In the latter case, the first term 
exponent has to be assumed. 

Should it be required to develop the equation for a curve passing 
through two points, in addition to the origin, two possible methods are avail- 
able. If the ordinate of the second point is larger than that of the first, 
and the slope is of no importance, a single term equation can be found, as 
shown in Fig. 13 and Eqs. (49), (50), and (51). If, however, a certain slope 
is specified at one or both points, or if the ordinate of the second point is 
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less than that of the first point, then the '"trial and error" method outlined 
in Section 2.60 must be used to find the desired equation. 

2.50 - Conversion from Two to Three or More Term Equations ;- In the previous 
sections, the development of two term equations was demonstrated. However, it 
may sometimes be found that, after a two term equation has been developed, and 
the ordinates computed and plotted, the resulting curve may not quite meet the 
requirements, and some slight modification is needed. This may be done by 
changing the exponents of one or both of the two terms, but such a method can 
involve several trial and error assumptions before the desired curve is ob- 
tained. 

A simpler method consists of adding a third term, the value of which 
can usually be found quickly. This method should not be confused with that of 
originally developing a three term equation for the curve, as will be outlined 
later on. Instead, the third term is added solely as a means of modifying a 
fundamental two term equation. 

Assume that the finally desired curve has already been drawn, in 
preliminary form - either sketched or faired with splines or french curves, 
or known from some previous design. The curve developed by the two term power 
curve equation, (which we have assumed needs modification), can then be placed 
over this preliminary curve, the points of origin made to coincide, and the X 
axis raised or lowered until the two term curve most closely coincides with the 
desired curve. The X axis of the two term curve will then form an angle with 
the X axis of the other. The tangent of this angle then becomes the constant of 
a third term in (x), the exponent of which is unity.- If the X axis of the two 
term curve has to be lowered the third term will be negative, and positive if 
the X axis is raised. 
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Consider, for example, the curve represented by Eq. (64). Let It be 
assumed that, in order to make this curve coincide more closely with a desired 
curve, it is necessary to raise the X axis by 0.5 inches at the Station (x) » 
100 inches. The constant of the third term will then be + 0.005, or, in log 
form, N(0. 69697-2). Hence, Eq. (64) would be modified to: 

y - N(0.66636)(x) * 4 + N(0.69897-2) (x) - N(0.28565-5)(x) 2 ' 625 (67) 

Changing the equation in this fashion, instead of initially developing 
a three or more term equation, will sometimes save time and computations. How- 
ever, when using this method, the «chord n , or distance from the origin to the 
point of intersection with the axis, will increase depending upon the amount 
the axis was inclined. To find the length of the new chord it will be necessary 
to compute, and plot, a few stations near where the point of intersection with 
the axis is expected. If this new chord Is not materially different from the 
original one, but it is still necessary to maintain exactly the original chord, 
then the last term should be changed, or an additional term added. 

2.60 - Development of Multi-Term Power Curve Equations ;- In Section 2.30, it 
was demonstrated how a single term power curve equation could be developed for 
a curve starting at the origin and passing through two given points, or through 
one point at which the slope was known. (See Figs. 13 and 14). However, when 
using the single term form, the curve could not be made to become parallel to, 
or intersect, the axis at any time, but instead continued ever further away 
from both X and T axes as (x) was increased. 

To make the curve of a single term equation become parallel to the X 
axis, and eventually intersect it, a second term can be derived and subtracted 
from the first, making a two term power curve equation as discussed in Sections 
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2.40 through 2.43. The two term power curve equation can only be developed 
when certain conditions are specified, however, and hence is somewhat limited 
in its use, although not nearly as much as the single term equation. 

Therefore, in order to derive a curve which can meet more specified 
conditions, or can be more generally used, a multi-term equation consisting 
of three or more terms must be employed. The development of such an equation 
may require several trial and error assumptions, but the final result will be 
a smooth continuous curve, with uniformly varying curvature, and passing 
through the desired points. The number of trial and error processes required 
to develop the equation will also be found to be substantially reduced after 
a little experience has been gained in the use of the method. 

To illustrate the method used in this process, consider Fig. 19, where 
it is assumed that a multi-term equation has to be found for a curve passing 
through points (P) and (Q), and some other point which need not be yet given. 
It will also be assumed that the required slope at point (Q) is given. The 
process will only be shown here for the first two points, as a numerical example 
later on will illustrate the entire development of a multi-term equation. 




FIGURE 19 
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A single term equation can easily be found for the solid curve in 
Tig. 19 passing through point (P), since, by assuming the exponent, (n), the 
constant (k) can be found from Eq. (53) and the equation written in the form: 

7 - kW D 
If the required ordinate and abscissa of point (Q) are B» and A*, 
respectively, substitution of this value of A 1 in the equation for the curve 
passing through point (P) will give the ordinate of this curve at Station R, 
or at the same abscissa as point (Q). Let this ordinate be called (E). The 
ordinate desired at point (Q) is B\ however, and therefore the second term of 
the equation will have to reduce the ordinate (E) by the amount (E - B'), or 

<C). 

Similarly, since the required slope at point (Q) is given, we can 
substitute the value of A 1 in the derivative of the first term and find the 
slope the solid curve will have at Station R. If this slope at point (Q) which 
is required is called (t'), and the slope of the solid curve at Station R is 
designated (t), the second term of the equation will have to change the slope 
by an amount (t - t 1 ), which we shall call (D 1 ). 

"Writing these equations: 

y. » k^A 1 )* 1 - C - B* ........(68) 

or, since E - ^(A 1 ^ 1 ? we get: 

C' -E-B' — ....-.. — (69) 

and: 

(dy/dx) = n^A'A- 1 - D« = f <™> 

/ t N D-, -1 - ' 

or, since t = n-jk^A ) J- : . 

D ' =t-t' ............(71) 
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The second term can then be considered as the equation of a curve 
which, at an abscissa of (A 1 ), will have an ordinate of (C 1 ) and a slope of 
(D 1 ). Hence, from Eqs. (55) and (54), the required constant and exponent, 
respectively, of the second term are: 

k 2 = C'AA 1 ) 112 (72) 

and: 

n 2 = (aV)/(C') ....(73) 

The equation has now been developed to the extent: 
y = k x (x) n l - k 2 (x) n 2 

In the above equation, (i^), (kg), and (n g ) are known, and (n-^ is 
also known since it was assumed at first from a consideration of the shape of 
curve desired at the origin. (See Figs. 15 and 16). 

In evaluating the equation at this point, it will usually be found 
that the ordinate at point (P) has been reduced by the second term, if this is 
negative, or increased if the second term is positive. Should such be the case, 
then a re-development of the equation becomes necessary with a correction to the 
value of the ordinate used in determining the constant and exponent of the 

first term. 

For the third term of the equation, the procedure to be used is the 
same a3 just outlined for finding the second term, providing the slope and 
coordinates of the third point are known. If one of these is unknown, it must 
be assumed in order to permit the slope and ordinate corrections to be found. 
Here again, after the first value of the third term has been determined, the 
entire equation should be evaluated to ascertain the effect of this third term 
on the ordinates. If the effect is too large, the entire process has to be 
repeated with corrections made to the first two terms, and a new third term 
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found. 

It should be remembered that by employing these trial and error 

corrections, and developing new equations to finally make the curve fit the 

specified points, the position of the maximum ordinate might shift somewhat. 

In many oases, this shift will be so slight as to be of no importance. If the 

position of the max!*™ ordinate must remain fixed, however, then perhaps 

another exponent can be assumed for the first term, or the abscissa and/or 

the ordinate of the first or second potot can be changed somewhat without 

affecting to any degree the overall requirements for the final complete curve. 

So far it has been assumed in most cases that negative terms are used 

after the first term. In many cases, though, positive terms are needed as for 

c^ple when reverse curvature is desired. The signs of the t*rms are usually. 

determined automatically in the development, however, for if (0 turns out to 

be negative when a negative term is assumed, the term should be positive, and 

vice versa. In addition, the signs can almost always be determined by visual 

inspection. 

m order to obtain a smooth continuous curve to pass through a given 
set of points, every attempt should be made to obtain exponents successively 
increasing in magnitude, such as 0.5 for the first term, from 1.0 to 9.0 for 

developed for a given set of points does not have exponents increasing in 
such an order, it indicates that the specified ordlnates of one or mere of 
the points are too large or small to fit on a smooth curve. Hence, the ordinate 
or equates in question should be modified slightiy and a new equation develop- 
ed. It will ataost always be found that such modifications can be made readily, 
as the ordinate of each point is usually specified as a maximum or minimum, 
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and therefore the correction can be made in the opposite direction. 

It will be noted that in the discussion on Fig. 19, the slope at 
point (Q) was assumed to be known. If the ordinate at point (Q) Is greater 
than at point (P), however, and the slope is not specified at point (Q), a 
single term can be used for the curve passing through these two points. In 
this case, the exponent would be found by Eq. (51), and the constant determined 
by substituting this exponent back in an equation similar to Eq. (49) or (50). | 

Thus, by developing another term for the last point on the curve, and combin- 
ing this term with the single term for the first two points, a two term equation | 
would be obtained for the entire curve, instead of the three term equation < 

obtained in the example discussed. § 

' f 

The reader should carefully note the above statement that the use of I 

a single term for the first two points only applies when the ordinate of the 1 

second point is larger than that of the first. This is because a single term 

equation, as noted previously, can not be made to define a curve which at some 

point becomes parallel to the X axis and then reverses its slope. Hence, should 

the ordinate of the second point be smaller than that of the first point, and 

the slope of the second point is unknown, then this slope must be assumed and 

two terms found, as before, for the first two points. (Eqs. (68) through (73)). 

An illustrative example is given below to show how the equation can be 

developed for a curve passing through three points. Suppose that the curve is | 

to pass through the three points, (P), (Q), and (R) as indicated on Fig. 20, 

and that the slope at point (Q) is to be zero, and the slope at point (R) is | 

to be -0.22. f 

An exponent of 0.5 will be assumed for the first term. Hence, for % 

point (P): 



FIGURE 20 
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A - 20 
B = 8.5 

Since we know that the second term will have to be negative in order 

to make the curve intersect the X axis, it will reduce the value of (B) at 

point (P). Therefore, assuming for the first term that (B) =9.0, (k) can he 

found from Eq. (53): 

k = B/(A) D 

=9.0/(20)°* 5 

or, in logs: 

k = log(9) - (0.5)log(20) 

* N(0.95424 - 0.65052) 
» N(0.50572) 
Therefore the first term is: 

y = N(0.30372)(x) * 5 ....... (74) 

The value of (7) when (x) = (A*) = 40 is: 
7 =N(0. 30372) (40) 

= N(0. 30372 + 0.5 log 40) 
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= N(0. 30372 ■*- 0.80103) 
= N(l. 10475) 
= 12.728 
The correction which has to be applied to the ordinate by the second 
term is then, from Eq. (69): 

C* = E - B' 

= (12.728 - 10) 
. - 2.728 
The slope of the curve defined by Eq. (74) at an abscissa (x) = (a') = 
40 is found from the first derivative: 

(dy/dx) = (0.5)N(0. 30372) (40)"°* 5 

= N(log 0.5 •»- 0.30572 - 0.5 log 40) 
= N((0. 69897-1) + 0.30372 - 0.80103) 
= N(0.20166-l) 
= 0.15910 
The slope correction to be applied by the second term is found from 
Eq. (71), where, in this case, (t ), the required slope at point (Q), is zero: 

D = t - t 

= (0.15910 - 0) 
= 0.15910 
The exponent of the second term is now found from Eq. (73): 

; n 2 =(AV)/c' 

= (40)(0.15910)/(2.728) 
= 2.3328 
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and the constant is found from Eq. (72): 

k 2 = c7(A*) n 2 

= 2.728/(40) 2 - 3328 

or, in logs: 

kg - N(log 2.728 - 2.3328 log 40) 






= N(0,45584 - 5.73733) ^# 









s N(0.69851-4) } : r'.H 

Since C* and D* are positive, the second term roust be negative, and the 

*^L.jls*r e 



equation so far becomes: ' ": 
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y = N(0.30372)(x) - 5 - N(0.69851-4)(x) 2 ' 5528 (75) 

The third point has the ordinate A* = 100, B n - 0, and the slope is ' " y 
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-0.22. To determine the ordinate correction to be supplied by the third term, ,"'; 
the two term equation as given by Eq. (75) must first be evaluated for (x) = 
A = 100: ■■ ■ 'Afrjn* 

y = N(0. 30372) (100) * 5 - N(0. 69851-4) (100) 2 - 3328 

s N(0.30372 + 0.5 log 100) - N( 0.69851-4) (2.5328 log 100) 

« N(1.50572) - N(1.36417) 

s 20.124 - 23.129 

- - 3.005 .....(75a) 

■:■/.. 

■■■,. ; . " - _ *■ -.;* u r 

The minus sign indicates that the curve defined by Eq. (75) would lie 
3.005 units below the X axis at (x) = A* « 100, and therefore the third term 
should be positive in order to correct this condition. Hence, from Eq. (69), 

■-'■■-.■ . v- s- *«&: 

■ '.■'.-■ ■* * -w x & 

the correction to be applied to the ordinate by the third term is: - \.'%\7 

C = E - B v ; '^8>. 

, (-3.005) - < r ^&P 

= -3.005 ' (75b) r ~ 
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The slope correction to be applied by the third term is found by first 
determining the slope Eq. (75) would have at (x) = A » 100. Therefore, differ- 
entiating Eq. (75), we haves 

(dy/dx) - (0. 5 )N (0.30372) (x)" * 5 - (2.3528)N(0.69851-4)(x) 1 - 3328 

For (x) » A » 100, this equation has the values 

(dy/dx) = N[log(0.5) + (0.50372) - 0.5 log (100)] - N|log (2.3528) + 
(0.69851-4) + 1.3328 log (100)] 
= N(0.00269-l) - N(3. 73205-4) 
= 0.10062 - 0.53958 

= - 0.43896 (75c) 

Here again the sign is negative, showing that the third term must be 
positive to raise the slope from -0.43896 to the specified -0.22. Hence, from 
Eq. (7l)r 

D - t - t 

= (-0.43896) - (-0.22) 

= - 0.21896 (75d) 

The exponent of the third term is now found from Eq. (73): 

n s = A D /C 

= (100) (0.21896 )/(5.005) 
= 7.28650 

The constant is found from Eq. (72): 

k 3 = c7(aY 3 

= (3.005)/(100) 7 * 28650 



Iff 



= N [(log 3.005) - (7.28650 log 100)] 
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- N(0. 47784 - 14.573) 

= N(0. 90484-15) 
Note that in determining the values of (n s ) and (k 3 ) above, the 
negative signs for c" and d" were not used, since they are only used to determine 
the sign of the constant (kg), which has been shown to be required to be 

positive. 

The equation of the curve is now found by adding the third term to 

Eq. (75): 
7 = N(0.30572)(x)°- 5 - N(0.69851-4)(x) 2 - 5528 + 1 ( 0. 90484-15 )(x) 7 - 2865 ° ...(76) 

in computing the ordinates from Eq. (76), it. will be found that at 
point ? the ordinate is 8.459 instead of the specified 8.5. The difference is 
0.041, which, in most cases, is negligible. Should a closer result be needed, 
however, the equation should be re-developed starting with an ordinate at 
point P of perhaps 9.1 instead of the value of 9.0 assumed in this example. 

2.61 - ^termination of Terms w hsn Slopes are Unknown: - In the 
illustrative example given in Section 2.60, it was assumed that the slopes at 
points Q and R, (Fig. 20), were known. Very often, however, the slopes are not 
known, and the equation must be developed in another manner. 

As mentioned previously, if the ordinate at point Q is larger than that 
at point P, and the slope is unknown, a single term can be found to fit these 
two points by the use of Eqs. (51) and (49) or (50). To then find the second 
term, (for the third point), the slope at the third point must be assumed, and 
the procedure followed as given in the illustrative example in Section 2.60. 

Another method by which the equation may be developed for three points, 
when all the slopes are unknown, is to first find a single term for the first 
point by either assuming a slope at this point or the exponent of the first 
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term, Eqs. (53) and (54). Then this single term equation can be evaluated up 

to the abscissa of the third point, and a second term found to take into account 

the last two points, as demonstrated below. 

Consider the illustrative example shown by Fig 20 in Section 2.60. Let 
it be assumed that, in addition to passing through points P and Q, the curve 
must also pass through the points (x) ■ 90, (y) » 2.5 and (x) =■ 100, (y) =0, 
(point R). Here we have a case where four points must be satisfied, but it will 
be noted that the slopes at the last two points are not specified. This example 
is then indicative of the condition mentioned above where a term is needed to 
take into account the last two points on the required curve. 

Eq» (75) gives the equation as developed for the first two points, P and 

Qr 

y = N(0.30372)(x)°* 5 - N ( 0. 69851-4) (x) 2 " 532 ! (75) 

This equation must now be evaluated at the abscissae at the last two 
points, (x) = 90 and (x) - 100. Hence, for (x) = 90, the value of (y) for Eq. 
(75) is: 

y -N(0. 30372 + 0.51og 90 ) - N(0.69851-4 + 2.5328 log 90 ) 
= N(0.50372 + 0.97712 ) - N(0.69851-4 -j- 4.55891) 
= 19.091 - 18.089 
. 1.002 

At station (x) = 100, the (y) for Eq. (75) is: 

y a N(0.30372 + 0.5 log 100) - N(0.69851-4 + 2.5328 log 100 ) 
= N(1.30372) - N(5.36417-4) 
= 20.124 - 23.129 
= - 3.005 
The ordinate correction at station (x) = 90 is, from Eq. (69): 
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C = E - B 



= 1.002 - 2.50 
= - 1.498 
and, at station (x) = 100, where (y) = E = 0: 



C"' - E - B W 



s (-5.005) - . 
■ ss - 5.005 
Both of these corrections have negative signs, showing that the third 
t«- «t be posits. Hence, at station M - 90, the reared constant of the 
third term is, from Eq. (72): 



,n J, ,H\n? 



k 3 = c7(A") n 3 



(l.498)/(90) n 3 



(77) 



and, at station (x) = 100: 



k 5 = c7(A") n 3 

= 5.005/(l00) n 3 -(78) 

The value of (n 5 ) required to satisfy hoth B*u (77) and (78) can no* 

be found by equating the two: 

(1.498)/(90) n 5 = (5.005)/(100) n 3 

(100/90) n 3 = (5.005/1.498) 

(l.Ul) n 3 = 2.007 

n« = log(2.007)/log(l.Ul) 
5 

s 6.61695 
The valne of <*fc) is no. fovmd by snbstitnting back i» either E,. (77) 
or (78), as follows: 



k* = 3.005/(100) 
o 



6.61695 
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s N[log(3.005) - 6.61695 log(lOO)] 
s N(0.47785 - 15.23390) 
= N(0.24595-15) 
The third term is then: 

y = N(0.24395-13)( X ) 6 * 61695 
and the final equation is Eq. (75) plus this third term, ort 

y - N(0.30572)(x) - 5 - N(0.69851-4)(x) 2 ; 3528 + N(0.24395-15)(x) 6 * 61695 .....(79) 

The above shows the manner in •which equations can be developed without 
the slopes being known, or the method by which two points can be used to permit 
the development of a single term in the equation. 

By employing the principles evolved in this and the proceeding section, 
multi— term power curve equations can be developed for a curve passing through 
many points. Obviously, the greater the number of points specified through 
which the curve must pass, the more complicated the equation becomes, involving 
more and more computations. Therefore, it is advisable to keep to a minimum the 
number of points required to lie on the curve. It will usually be found that this 
CJkU be done in many instances, for although it may be initially specified that 
a ll the given points should lie on the curve, many of these points will be of 
secondary importance, or will be maximum or minimum values where slight variat- 
ions can be permitted. 

2.70 - Correction of Equations for Change in Thickness Ratio s- The ratio of the 
external dimensions of thickness and length of a body are of considerable use 
and interest where any kind of fluid flow is involved. This ratio is known as 
"thickness ratio", and is defined as the maximum ordinate of the body, or curve, 
divided by the length. This is usually expressed in percent of the length. Thus, 
for example, if the maximum ordinate is 12.5 units and the length 100 units, the 
thickness ratio would be 12.5/100, or 0.125, or 12. 5%. 
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It Is very often desired, after the equation has been developed for a 
desired curve, to change the thickness ratio for one reason or another. This 
can be most easily done, -without resorting to the necessity for developing 
another equation from the beginning. 

Consider, for example, Eq. (64) in Section 2.42, where the maximum 
ordinate was specified as 26 inches and the length, or chord, as 262 inches: 

y - N(0.66656)(x)°* 4 - N(0.28565-5)(x) 2 ' 625 ........... .....(64) 

Suppose that it is desired to change the maximum ordinate to 25 inches, 
and a new equation Is therefore needed. It has been previously noted that the 
constants of an equation control the thickness ratio, and that the exponents 
have nothing to do with it. Therefore, if the thickness ratio is varied, only 
the constants have to be changed, and this occurs in the same proportion. 

Hence, to change from a thickness of 26 inches to 25 inches in the above 
example, without changing the chord, it is only necessary to multiply the 
constants by the ratio 25/26, or, using logs, the correction to be added to each 
constant becomes: 

= log(25) - log(26) 
■ ' > - = (2.39794-1) - (1.41497) 

~ 0.98297-1 . ' „ 

Adding this amount to the log of each constant in Eq. (64), the 
corrected equation becomes: 

y = N [(0.66656) + (0.98297-1)] (x)°* 4 - N [(0.28565-5) + (0.98297-1)] (x) 2 * 625 
= N(0.64933)(x)°- 4 - N(0.26862-5)(x) 2 - 625 (80) 

Should it be desired to thicken a curve, the same procedure is used, 
the multiplier of the constants then being greater than unity. In the above 
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example, a two term equation was used for convenience, but the method is 
applicable to equations with any number of terms, provided, of course, that 
yll the constants are corrected in the same ratio. 

To change the thickness ratio by lengthening or shortening the chord, 
while leaving the value of the maximum ordinate unchanged, the procedure becomes 
slightly more elaborate. In this case, all the (x) ? s have to be multiplied by 
a number equal to the ratio of the chord on hand to the desired chord. Since 
these (x)'s are raised to a power, the multiplier must also be raised to the 
same power. The multiplier, however, being a constant, can be separated and 
raised to the desired power and then can be included in the constant of the 
term. 

Considering Eq. (64) again, assume that it is desired to change the 
chord from 262 inches to 270 inches, leaving the m a xim um ordinate at 26 inches 
as in the original curve. Then each (x) must be multiplied by the inverse ratio 
of the change in chords, or 262/270. Using logs, the multiplier becomes: 

= log (262) - log (270) 

- (3.41830-1) - (2.43136) 

= 0.98694-1 

Substituting this value in Eq. (64), we get: 

y = N [(0.66636) + 0.4(0.98694-1)] (x)°- 4 - H [(0.28565-5) + 



2.625(0.98694-1)] (x) 2. 



625 



To demonstrate the method of multiplying logarithms by fractional 
constants, the above equation will be evaluated in detail: 

0.4(0.98694-1) = (0.39478-0.4) 



^ .(jaaaB^swmsaMwww&jii 



65. 

The anti-log of the number on the right hand side of the above 
equation can not be found unless the negative number is a whole digit, such 
as 1, 2, 5, etc. Hence, 0.6 can be added to both numbers in order to make the 
negative number equal to -1, or: 

= (0.39478 + 0.6) - (0.4 ■»- 0.6) 
= 0.99478-1 
This logarithm is now in a form to be added to, or subtracted from, 
other logarithms to permit the final anti-log to be found in log tables. 
In a similar manner: 

2.625(0.98694-1) = (2.24070-2.625) 
Adding 0.575 to both numbers to make the negative number a whole digit: 
= (2.24070 + 0.375) - (2.625 + 0.575) 
= 2.61570-3 
These values must now be added to the respective constants in Eq. (64), 
as shown below: 

The constant for (x) * becomes: , 

= N [(0.66636) + (0.99478-l>] = N(0, 66114) 

, v2.625 , 
and the constant for (x) becomes: 

= N [(0.28565-5) + (0.61570-1)]= N(0.90135-6) 
Eq. (64), corrected to a chord length of 270 inches and a maximum 
ordinate, as originally, of 26 inches, then is: 

y = N(0.66114)(x)°* 4 - N(0.90135-6)(x) * (81) 

Note that, by this method of including the multiplier in the constants, 
the exponents are unchanged, thus keeping to a minimum the number of computations 
required to compute the ordinates of the new curve. 
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CHAPTER III 
CURVATURE AND POINTS OF INFLECTION 

3.00 - Points of Inflection on a Curve ;- It has been previously stated that 
the curves defined by one term power equations start at the origin, turn 
swiftly in one direction or the other, (depending on the value of the exponent), 
and then continue their course with the ordinate continuing to increase as the 
abscissa is increased. The direction of curvature is always the same, or, in 
other words, the curve continues to bend more and more in the same, direction 
for as far as it is extended. Such a curve is known as a B continuous* curve, 
for, although it will never bend enough to become exactly parallel to either 
axis, it will continue to bend infinitesimal amounts as it approaches infinity. 

When a negative term is added to a single term equation with an exponent 
less than unity, the tendency is to make the curve bend around some more. If 
the exponent of the second term is larger than that of the first term, the curve 
will eventually become parallel to the axis and then intersect it later on. 
This curve also does not change it3 direction of curvature, and therefore is 
continuous. 

An easy way to understand the meaning of a "continuous" curve is to 
visualize segments of the curve being drawn by a compass. As long as the segments 
can be drawn with the compass center on the same side of the curve, the curve 
Is continuous. When the compass has to be shifted to the other side of the 
curve in order to draw the next segment, the curve becomes "discontinuous*, or 
reverses. 

A reverse curve is usually obtained from a power curve equation, with 
a first term exponent of less than unity, when a positive term is added, provid- 
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_* that the exponent of this positive term is large enongh. Ohvioualr, if 
the positive term has a small exponent, or constant, and the previous negative 
term, or t», has a large exponent, the effect of the positive term will 
Just he to slights decrease the rate of curvature without affecting the contin- 
uity of the curve. Hence, it can not generally he said that all positive terms 
,U1 cause a curve to reverse, hut instead the mathematical equations of the 

curTe must be examined to determine if such a condition crista. 

It ^ also happen that a negative term rill cause a curve to reverse. 

For example, referring to Fig. *. Curve represents a power curve equation 

of two terns, the first heing positive rith an exponent larger than unity, with 

the second term being negative. 
Y 




y FIGURE 21 

Here the negative term causes the curve to reverse, since the first 
term exponent was larger than unity. Curve A represents a continuous curve 
„hich could he formed from a two term equation rith a first term exponent less 
than urity and a negative second term. If a third term should he added to this 
two term equation, end the term he positive rith a ^rge exponent, Cerve B might 
result, with a reverse in it. 



68. 

Many multi-term equations may have reverses in their curves so small 
as to be hardly noticeable by eye, and therefore it is necessary to mathematic- 
ally examine every equation developed to determine if any reverse in curvature 
is present. Reversals in curvature can be very detrimental to preserving a 
smooth fluid flow over a surface, and a breakdown of a smooth flow means an 
unnecessary drag, or resistance, for the object. Therefore, in streamlining 
work where fluid flow exists, and the resistance of the object to this flow 
must be kept to a minimum, reversals in curvature should be eliminated if at 
All possible, except sometimes near the end of the chord. 

To determine if any reverse curvature exists, the equation of the curve 
must be searched for "points of inflection", a "point of inflection" being 
defined as the point where the curvature changes its sign, or reverses. It will 
be remembered that the first derivative of an equation gives the equation of 
the slope of the curve at any point. 7/hen this equation is set equal to zero, 
and the value, or values, of the independent variable solved for, these values 
represent the points at which the slope is zero. The second derivative of the 
original equation, or the derivative of the first derivative equation, is a 
function of the curvature, and, when set equal to zero, gives the value, or 
values, of the independent variable at which the curvature is- zero. If the 
curvature changes sign after this point, this is a point of inflection as the 
curvature has reversed. 

The easiest means of finding whether such points exist is to plot the 
second derivative equation. Where the curve intersects the X axis the curvature 
obviously changes sign, and hence a point of inflection exists. Consider, for 
example, Fig. 20 and Eq. (76), where the specified slope at the end of the chord, 
(point R), was -0.22. The first derivative of Eq. (76) is: 
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v l. 3328 



*/*= (0.5)..(0.3037 S )W- - 5 - (*.5W8)««M9BBUOM 

,,., .6. 28650 
(7.2865)N(0.90484-15) (x) 

n!/ %-°' 5 nhnp 2 3328) + (0.69851-4) (x) 
= N[(log 0.5) + (0.30372)J(x) - H [(log 2.53^> + V j 

+ N [(log 7.2865) + (0.90484-15)] (x) 



1.3528 



.6.28650 



= N(0 



The 



.00 2 69)W-°- 5 - NtO-OaS^itx) 1 - 3323 + «<0.W™M*>M 



6.2865 



(82) 



second derivative is obtained in the same warmer: 
d 2 y /dx 2 = -(0.5)N(0.00269)(x)- 1 - 5 - (l.35 2 8)N(0.06659-5)(x) 

L 5.2865 



0.3328 



(6.2865)N(0.76736-14)(x)* 



L 5.2865 



-N(0 



been evaluated for varying values of 



Eqs. (76), (82), and (83) have 

Table IX and plotted on Fig. 22. 



(x), and the results are given in 



TABLE IX 
Evaluation of Eqs. (76), (82), and (83) 
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7 (Eq. 76) 




FIGURE 22 



Fig. 22 shows that the curve of the second derivative intersects the 
X axis at approximately (x) = 88 inches, indicating a point of inflection, or 
a point where the curve reverses. As a matter of interest, the point at which 
the slope of the original curve is zero was specified at (x) =» 40 inches, (see 
Fig. 20), and this is shown on Fig. 22 where the curve of the first derivative 
intersects the axis. This is the maximum ordinate position. 

Let us now assume that it is desired to change the equation so that 
the curve does not reverse at any point. This can be done by a trial and error 
method, either by assuming the slope at the end point or by assuming the ex- 
ponent of the last term. As a first approximation, therefore, let the slope at 
the end point, (x) * 100, be changed from the original value of -0.16 to -0.22. 
Only the third term need be changed, and the required corrections to be applied 
by the third term to the ordinate and slope at (x) = 100 can be found as 
originally done in developing Eq. (76). 

The ordinate correction remains unchanged from Eq. (75b), since no 
change has been made. Hence: 
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C» =-5.005 ( 75b ) 

The slope which the first two terms of the equation would have at 
(x) = 100 remains -0.43896, from Eq. (75c), but the slope correction is now, 
from Eqs. (71) and (75c): 

d" = -0.43896 - (-0.16) 
. » -0.27896 
The new exponent is then, from Eq. (73): 

n 5 = A D /C 

= (100)(0.27896)/3.005 
s. 9.28320 
and the new constant is, from Eq. (72): 

k5 = c7(AY3 

- 5.005/(l00) 9 ' 28520 
= N(0.91144-19) 
Eq. (76) is now, modified for a slope of -0.16 at (x) - 100 inches: 

y . N(0.50372)(x)°- 5 - N(0.69851-4)(x) 2 - 3528 + N(0.91144-19)(x) 9 - 2832 ? (84) 

Checking the effect that this new third term will have on the maximum 
ordinate at (x) = 40 inches, we substitute this value of (x) in Eq. (84) and 
find that the maximum ordinate is 10.0006 inches instead of the originally 
specified 10.00 inches. This effect is negligible and need not be further 

considered. 

To find the first and second derivatives of Eq. (84), only the third 
term need be investigated, as the first two terms will be unchanged from Eqs. 
(82) and (83), respectively. The derivation will not be carried out again, as 
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the procedure is the saae as shown in deriving Eqs. (82) and (85). Therefore: 
dy/dx = N(0.00269)(x)°* 5 - N(0.06639-3)(x) : 
and: 



.1.3328 



8.2832 



+N(0.87914-18)(x) ...(85) 



-1.5 



dv/dx =-N(0. 70166-1) (x) -N(0.19116-3)(x) 



0.3328 



+ N(0.79734-17)(x) 7 '^.( 86) 



Eqs. (84), (85), and (86) have been evaluated in Table X below, and 
plotted on Fig. 23. 



TABLE X 
Evaluation of Eqs. (84), (85), and (86) 



X 


7 


dy/dx 


d^y/dx 2 


1 


2.012 


1.0050 


-0.50466 


2 


2.841 


0.7086 


-0.17983 


5 


4.479 


0.4400 


-0.04765 


10 


6.257 


0.2931 


-0.01925 


30 


9.629 


0.0753 


-0.00788 


50 


9.644 


-0.0710 


-0.00699 


70 


6.882 


-0.2006 


-0.00554 


90 


2.151 


-0.2463 


+0.00319 


100 


0.000 -0.1600 


+0.01542 



Mg. 25 shows that, by changing the slope of the euro at station 
« = 100 inches fro., -0.22 t „ -0.16, the point at which a reverse occurs in 
the curve of the equation has shifted further forward than when the slope was 
-0.22 at the end point. This shift could have been seen beforehand, as obviously 
a decrease in the slope at the end point would -ove the position of reversal 
forward, but the example has been worked out to show the trial and error .anner 
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in which the problem must be solved. 

To eliminate the reversal in the curve, therefore, the slope at station 
(x) = 100 inches must be increased in a positive direction. For a second assumpt- 
ion In this case, then, the exponent of the last term will be assumed instead of 
the slope in order to demonstrate how this method may be used. ,,f -; : 




Point of Inflection 
dV/dx 2 (Eq. 86) > 

FIGURE 23 
the value of (d^y/dx 2 ) for the first two terms of Eq. (85), at station 

(x) = 100. inches, is: 

d^y/dx 2 = -N(0.70166-4) - N(0.85682-3) 
: = -0.007692 

In order for (d 2 y/dx 2 ) for the complete equation to be zero at station 
(x) = 100, the (d^y/dx 2 ) of the third term must be -0.007692. Assume now that 
the third term exponent is 5.435. Then, since the ordinate correction required 
of the third term is still -3.005, (Eq. 75b), the value of k 3 is found from ^ 

Eq. (72): 

- n // . n \n* 
k 3 - C /(A ) 3 

- 3. 005/ (100) 5 ' 435 

= N(0.60784-ll) 
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and the third term is: 



Hence: 



7 » N(0.60784-ll)(x) 5 * 455 (87) 



dy/dx - N(0.54304-10)(x) * (88) 



and: 

dVdx^ N(0.98993-10)(x) 5 ' 455 (89) 

Evaluating Eq. (89) for (x) = 100, we get: 

d^y/dx 2 = N(0. 85993-10) 
= 0.007243 

This does not yet quite check the required value of 0.007692, but for 
this example will be taken as a close enough check. The third term of the equat- 
ion as given by Eq. (87) will, however, eliminate any reversal even beyond 
station (x) » 100 inches, and can be therefore considered satisfactory. The 
final equation is now: 

y = N(0.30372)(x)°* 5 - N(0.69851-4)(x) 2 * 5328 + N(0.60784-H)(x) 5 ' 455 (90) 

The first and second derivatives of Eq. (90) are: 

dy/dx = N(0.00269)(x)"°' 5 - N(0.06639-3)(x) 1 ' 3528 + N(0.54304-10)(x) 4 * 43 ?..(91) 

and: 

n r, 15 3328 3.435 

d'V/dx 4 = -N(0. 70166-1) (x) -N (0. 19116-5 )(x) ' -f-N(0.98993-10)(x) ...(92) 

Eqs. (90), (91), and (92) are eveluated below in Table XI, and plotted 
on Fig. 24. 

Fig. 24 clearly shows no point of inflection until station (x) = 100, 
and therefore no curve reversal. Eq. (90) then represents the equation of a 
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TABLE XI 
Evaluation of Eqs. (90), (91), and (92) 



X 


7 


dy/dx 


d^y/dx 2 


1 


2.012 


1.0050 


-0.50466 


2 


2.841 


0.7086 


-0.17985 


5 


4.479 


0.4400 


-0.04765 


10 


6.259 


0.2951 


-0.01925 


50 


9.655 


0.0761 


-0.00776 


50 


9.708 


-0.0645 


-0.00647 


70 


7.204 


-0.1815 


-0.00512 


90 


2.697 


-0.2605 


-0.00250 


100 


0.000 


-0.2757 


-0.00045 




100 



7 (Eq. 90) 

dy/dx (Eq. 91) 

dVdx 2 (Eq. 92) 



Point of Inflection (slightly 
past (x) s 100) 



FIGURE 24 
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continuous curve, with no reversals, passing through the points (x) = 20, (y) = 
8.5 and (x) = 40, (y) » 10, with the maximum ordinate at station (x) = 40. 
These were the conditions specified for Fig. 20 and Eq. (76), except that it 
was found necessary to change the slope at the end point, (station R in Fig. 
20), in order to eliminate a reversal in the curve. This results in the new 
third term given in Eq. (90). 

3.10 - Determining the Smoothness of a Curve ; - It has been shown in 
Section 5.00 how the curve of the second derivative is used to check for points 
of reversal in the basic equation. Such reversals, as has been pointed out 
before, are usually very detrimental to the preservation of an undisturbed 
fluid flow over the surface, and can be the cause of excessive resistance of 
the object when moving through a fluid. In aircraft work, where the fluid speeds 
are continually increasing as aircraft speeds increase, the drag caused by 
reversals in the contour lines of fuselages, nacelles, etc., may be many times 
magnified. 

It is usually not sufficient, however, to only examine a curve for 
points of reversal, as abrupt changes in slope, or curvature, of the contour 
lines may have the same effect as reversals in curvature. It has long been 
known that the pressure distribution over a body is a function of the degree 
of curvature of the body, and, in turn, that this pressure distribution governs 
the magnitude of the local velocity of the fluid moving past the body. Rapid 
increases or decreases in curvature of the body cause rapid accelerations, or 
decelerations, in the local fluid velocities, resulting in a breakdown of the 
smooth flow, with corresponding high drags. 

It is therefore necessary, particularly where high fluid velocities 
prevail, to insure that the curvature of the body changes as gradually as 
possible. This can be checked by the use of the same second derivative curve 
used in examining for points of reversal, for this curve really represents the 
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rate of change of slope of the basic curve, and 1. therefore essentially a 
Action of its cu^e. Hence, a plot of the (*W) — ""*» ^ 

«*+>i no Deaks. or sudden changes, or wavmess. 
a 3 tnooth continuous line, wxth no peaics, or 

4,3= «-f* thp axis, for to cross the 
Obviously this curve should also lie on one side of the arts, 

axis would signify a reversal in the basic curve. 

Fo r example, referring to the (*,/*«> — » ««• 25 > " ^ *" 
aeen that this curve foras a smooth line up until short* before it crosses 
ths> X axis at the point of reversal in the hasic curve. The curvature of the 
.aaic line can therefore he considered to chan g e very Xittle fro. ahout station 
M - 10 to station (x) - 70, indicate a good smooth hasic line. Ahead of 
station (x, - 10 the curvature changes very rapidly, .bich, in this iXlustrat- 

. nested, front section TO s desired for other reasons. Aft of station (x) - «, 
tt . change of curvature is ,uite large, and, even if no reversal .ere present 
in the hasic curve at this point, this rapid change in curvature night prove 

Hmnfa smooth undisturbed fluid flow, 
to be auite harmful to the preservation of a smootn 

*. - +h«+ the change of curvature aft of station 
Fig. 24, however, shows that the cnang- 

W _ 70 TO s substantially reduced .hen the reversal in the hasic curve TO s 
existed. Here- the (dV^> curve is a very smooth flat line, shoeing that, 
aft of station (x) = 10, as hefore, the curvature of the basic line hard* 
changes. This .sens that the fluid, in flowing around the contour line, -31 at 
no pl.ce encounter sudden Ureases or decreases in curvature, and therefore 
should not suddenly brea* d 0TO and cause undue drag increases. 

m Xines developed, therefore, should he carefully checked for the 
possibility of sudden curvature changes, particularly .here high fluin velocit- 
ies are involved. Tf such changes in curvature occur, a study of the first and 
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second derivative curves will show where corrections should be made to the 
power curve equation. Obviously, no set rules can be stated for applying these 
corrections, but a little experience in such matters will enable the engineer 
to accomplish it in a short time. 

5.20 - Curvature and Radius of Curvature :- In Section 5.10 it was stated 
that the curve of the second derivative, (dry/dsr), can be used for checking 
the curvature of the basic line. This is quite true, but mathematically, the 
value of (d2y/dx2) does not give the actual curvature, but instead is a direct 
function of the curvature. 

Very often, however, it may be desired to know the actual curvature, 
or the radius of curvature. The latter is very useful in permitting the line 
to be reproduced, as it is sometimes easier to work with radii than with a 
table of ordinates. The expressions for curvature, K, and radius of curvature, 
R, are therefore given below. The derivations of these equations will not be 
outlined here, as they may be found in any good reference to the Calculus: 

K = (dVax 2 )/[l + (dy/dx) 2 ] 1 ' 5 (93) 

and: 

R * 1/K - [l + (dy/dx) 2 ] 1 - 5 /(d 2 y/dx 2 ) (94) 

As an illustrative example in the use of the expressions for curvature 
and radius of curvature, Eqs. (84) and (90) are evaluated below in Tables XII 
and HII, respectively. The calculations are self-explanatory, the values of 
(dy/dx) and (d2y/dx 2 ) for Eq. (84) being taken from Table X, and the correspond- 
ing values for Eq. (90) from Table XI. 

The values from Tables XII and XIII are plotted on Figs. 25 and 26, 
Fig. 25 showing the curvature, K, and Fig. 26 the radius of curvature, R. 
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10 



30 



50 



70 



90 



TABLE XII 
Curvature and Radius of Curvature for Eq. (84) 



100 



2.012 



2.841 



4.479 



6.259 



9.629 



9.644 



6.882 



2.131 



0.000 



dy/dx 



1-0050 



0.7086 



0.4400 



0.2931 



0.0753 



-0.0710 



-0.2006 



-0.2463 



-0.1600 



&y/dx* 



-0.50466 



-0.17983 



-0.04765 



-0.01925 



-0.00788 



-0.00699 



-0.00554 



+0.00319 



+0.01542 



K 



-0.17709 



-0.09769 



-0.03654 



-0.01701 



-0.00781 



-0.00693 



-0.00522 



+0.00233 



+0.01480 



R 



-5.64 



-10.22 



-27.40 



-58.70 



-128.00 



-144.30 



-191.50 



341.50 



67.50 



TABLE XIII 
Curvature and Radius of Curvature for Eq. (90) 



X 


7 


dy/dx 


dV/dx 2 


K 


R 


1 


2.012 


1.0050 


-0.50466 


-0.17709 


-5.64 


2 


2.841 


.. 0.7086 


-0.17983 


-0.09769 


-10.22 


5 


4.479 


0.4400 


-0.04765 


-0.03654 


-27.40 


10 


6.257 


0.2931 


-0.01925 


-0.01701 


-58.70 


30 


9.633 


0.0761 


-0.00776 


-0.00769 


-130.00 


50 


9.708 


-0.0643 


-0.00647 


-0.00642 


-155.90 


70 


7.204 


-0.1815 


-0.00512 


-0.00488 


-205.00 


90 


2.697 


-0.2605 


-0.00250 


-0.00233 


-430.00 


100 


0.000 


-0.2757 


-0.00045 


-0.00040 


-2500.00 



^'f f*? °^f ^■■r-^rr - ^-'-" "-■'"'--"' " 'r ^r- 1 ^^* yyrit r'ffiB^tt l i ffn'- : - ,-^Yr - " n — • -T r - " ^*^" *^-" a ^ j *- t ■^-j-^-^--^" -^^.:^- - ..~~~, 



80 i 



Basic Curve (Eq. 90) 



Basic Curve (Eq. 84) 



Curvature 
(Eq. 84) 



LOO I 




FIGURE 25 




Radius of 
Curvature 
(Eq. 84) 



*"'W- 



Radius of Curvature 
(Eq. 90) 



FIGURE 26 



"■""■■^ *r*xa»x:xx:4etngaif.z.... 



81. 



An inspection of Fig. 25 shows that the lines of curvature for 
Eqs. (84) and (90) are substantially the same as the (d ? j/dx Z ) curves of 
Figs. 23 and 24, respectively. This indicates why the (d^/dx 2 ) curve can 
be used in most cases for checking the smoothness of the basic curve. However, 
when the values of (dy/dx) are large, or vary considerably, it is advisable 
to use the actual values of curvature, and not (dSy/dx 2 ), in checking the basic 

curve. 

Fig. 26 clearly shows the reversal in curvature in Eq. (84), as it will 
be noted that the radius has shifted from one side of the X axis to the other. 

Further uses for the values of curvature and radius of curvature will 
be outlined in later sections. 
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CHAPTER IV 
FURTHER APPLICATIONS OF THE POWER CURVES 



The previous chapters have shown the methods whereby power curve 
equations can be developed to obtain various curves and shapes meeting certain 
requirements. In addition, several numerical examples have been given during 
th8 course of these developments illustrating how the equations are used in 
these cases. 

In this chapter, further examples will be given to show the application 
of the power curve equations. Quite obviously it is not possible to cover all 
of the conditions under which the equations can be used, but a thorough study 
of the examples already presented, and those to be presented in this chapter, 
will enable the reader to develop his own equations for any condition desired. 
4.00 - Combining a Power Curve with a Circle or an Ellipse ;- Power curves can 
be combined with a circle or an ellipse in a similar manner as an ellipse can 
be combined with a circle or another ellipse, a3 shown in Sections 2.21 and 2.22. 
For example, assume that it is required to combine a power curve with the circle 
shown in Fig. 27, and that the power curve is to have its origin at point and 
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be tangent to the circle at station C, 55 units from point 0. 

In order to determine the values of the exponent and constant for the 
power curve equation, the ordinates and slopes at the point, or points, through 
which the curve is to pass must first be found. Hence, since the curve is to be 
tangent to the circle at station C, the ordinate and slope at station C must be 

found. 

The ordinate may be found from either Eq. (2) or Eq. (7a), depending . 

upon where one- chooses to call the origin of the circle. Using Eq. (7a), the 
value of (y) at station C is: 

............(7a) 



./ 



2rx - x' 



= l/(2)(30)(25) - (25) 2 
• 29.580 
and, from Eq. (9), the slope at staion C is: 

dy/dx = (r - x)/(y) 

= (30 - 25)/29.580 
. 0.16905 
Developing xhe power curve ' equation for A - 55, B - 29.580, and t - 
0.16903, and using Eqs. (55) and (54): 

n - (55) (0.16903 )/29. 580 
= 0.S143 



(9) 



and: 



or: 



., .0.5143 
k = 29.580/(55) 

log(k) = log (29. 580) - 0.3143 log(55) 

= 1.47100 - 0.54700 

= 0.92400 

k = W(0.92400) 
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The required power curve equation is then: 
y = N(0. 92400) (x)°- 5143 

Now it may happen that the curve of the equation so developed is too 
blunt, or too sharp, or is unsatisfactory for some other reason. If such is the 
case, then a second point should be specified, say at a station 20 units from 
point 0, through which the curve must pass, and a two term equation found by 
the trial and error method outlined in Sections 2.60 and 2.61. If this is still 
not satisfactory, additional points should be specified, and a three, or four, 
or more, term equation developed by the same trial and error method. 

Consider, as a second example, the case shown in Fig. 28, where it is 
required to combine a power curve with an ellipse, the major axis of the ellipse 
being vertical. This problem can be solved by rotating the sketch so that the 

Y 




FIGURE 28 
major axis coincides with the X axis, thus being consistent with all the 
sketches illustrated so far, but actually it makes little difference to leave 
the sketch in the form of Fig. 28, as will be seen in the following development. 
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^posa tht « is m-x- *. Join a pee, C ^e .ith tna 1« 
elllp3 e at Station C. Here a g ain, tns slope ani tne atsoissa at station C 

have to be found first. 

fn _ ftf +hp pquation of the ellipse, 
Eq. (33) gives the power curve form of tne .qua 

. .,, * ^ +Vl the X axis and the origin at point B in 
the major axis being coincident wxth the X axxs 

«. -«r axis coincides with the Y axis, it is only necessary 
Fig. 28. "When the major axis coinciu«» 

to interchange (x) and (y), or: 



where, from Eqs. (32a) and (52b): 



x - (kjy - V S)2 



kl = (2b 2 )/a 



2w < 32a) 



8041 o ? .......(32b) 

k =(b 2 )/(a 2 ) 



2 

Hence: 



= (2)(30) 2 AO 



= 45.0 



and: 

k = (30) 2 /(40) 2 



- 0.5625 
it station C, (y) = 28, and therefore, from *. (33): 

x - /(45)(287^oT5625)(28) 2 
- 28.618 

4. + 4-.-nn C may be found, but first (x) and 
From Eq. (34), the slope at station C may tt. 

* a +>.p «res have been interchanged: 
( y ) must be interchanged since the axes hav. 
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dx/dy = (1^ - 2kgy)/(ac) * ( 34) 

_ (45) - (2) (0.5625) (28) 
(2) (28.618) 

= 0.2534 

The power curve equation may now be developed by again interchanging 
(x) and (y), and writing the equation in the form: 

x = k(y) n • -i 

If such were done, the exponent, (n), would be less than unity. However, 
the equation can just as well be developed in the familiar form, y= k(x) n , the ' 

exponent simply becoming larger than unity to give the vertical shape desired. 
(See Fig. 15). To do this, the slope, (dy/dx), is needed, and this is simply: 

dy/dx s l/(dx/dy) | 

- 1/0.2534 
= 3.9463 

Therefore, for t = 3.9465, B =* 72, and A = 28.618, Eqs. (53) and (54) < 

give: 

n - (28.618) (3. 9463)/(72) 
= 1.5686 

k - (72)/(28.618) 1 * 5686 

= NJlog (72) - 1.5686 log(28.618)] -J 

= N(0. 58275-1) ' : 

The required apation is then: | 

y = N(0.58275-l)(x) 1 * 5686 ? 



1 %, 
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In the combining of a power curve with an ellipse or a circle, or any 
other shape, to become tangent at any point, there is a limit to the extent to 
which this can be carried without encountering a reversal in the curves. This 
limit is that the origin of the power curve must lie within the straight line 
extensions of the tangent lines of the points of junction. 

4.10 - Supplementary Curve s t - It often happens that it is required to 
develop a curve to join two other curves, or that a given curve should be ex- 
tended in such a manner that, for example, one section of the extension should 
be a certain shape, such as an ellipse or circle, and the other section another 
shape, such as could only be defined by a power curve. For example, consider 
Fig. 29, in which it is assumed that the curve to the right of station 200 is 
to be extended to the left in such a manner that the section from station to 
station 40 is elliptical, with the ordinate at station 40 equal to 28. Between 
stations 40 and 200 another curve is to be developed which shall be tangent at 




FIGUEE 29 



station 40 to the ellipse and at station 200 to the original curve. The slope 
at station 200 is 0.09136, and the ordinate 54.876. 

It Is not possible, in this case, to develop the usual power curve 
equation, (with its origin at point 0), for the entire curve to station 200, 
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since it was specified that from (x) = to (x) = 40 the curve should be an 
ellipse. Neither is it possible to develop the usual form of power curve equat- 
ion with its origin at point P, since the slope at the origin of a power curve 
equation in the form y = k(x) n is infinite. In this case, the slope at point P 
of the curve from (x) = 40 to (x) = 200 must be the same as the slope of the 
ellipse at this point. 

Hence, the equation for the curve from (x) - 40 to (x) - 200 roust be 
developed in another manner. This can be done by first assuming the slope at 
point P. This slope can not be less than the slope formed by a line drawn 
through points P and Q, and must be smaller than the slope formed by a line 
through points P and 0. In the first case, the slope is: 

dy/dx = (54.876 - 28)/ (160) 
- 0.1695 
and, in the second case: 

dy/dx = 28/40 
= 0.70 
Therefore, let the slope at point P be assumed as 0.30. Now if a 
straight line be drawn thuough point P with this slope, and the origin taken 
at-point R, its equation is: 

y - 28.0 + 0.50(x) (95) 

and, at station 200, (or (x) = 160 from point R, the origin for this 
section of the curve), the value of (y) would be: 

y - 28.0 -»- 0.30(160) 
- 76.0 
This is larger than the required value of 54.876, and therefore a 
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correction must be applied in the manner demonstrated in Section 2.60. It will 
be noted that this development is exactly the same as outlined in Section 2.60, 
except that here we start by applying corrections to a straight line in order 
to curve it back to a specified ordinate, whereas in Section 2.60 the correct- 
ions were applied to an already curved line whose equation was in the form y » 

k(x) n . 

From Eq. (69), then, the correction to the ordinate at (x) - 160 is: 

C' = 76 - 54.876 
= 21,124 

To find the slope correction necessary, the first derivative of Eq.. (95) 
must be found and then evaluated at (x) = 160: 

dy/dx = + 0.30 
= 0.30 

The slope is obviously constant, since Eq. (95) represents a straight 
line, and the slope correction required at point Q is, fron Eq. (71): 

d' = 0.30 - 0.09136 
= 0.20864 

Since c' and D* are positive, the required correction term is negative. 
From Eqs. (72) and (73) the required constant and exponent, respectively, are 

now found: 

n = (160)(0.20864)/21.124 

= 1.5803 

and: 

N 1.5803 
k = (21.124)/(160) 
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= N [log(21.124) - 1.5803 log(l60)] 
» N (4. 33477-3 - 3.48317) 
= N(0. 85160-3) 

The equation for the curve from point P to point Q is then given by: 

1 5803 ' 
7-28.0+ 0.50(x) - N(0. 85160-3 )(x) * (96) 

It should be remembered that Eq. (96) assumes the origin at point R, 
and this must be allowed for in plotting the curve. 

For the elliptical part of the curve from point to point P, all that 
is known is one point on the ellipse, where (x) - 40, (7) -28, and (dy/dx) » 
0.30* It is therefore necessar7 to define the ellipse in terms of its semi-axes, 
(a) and (b). These values can be easil7 found, however, D7 the simultaneous 
solution of Eqs. (25) and (26): 



7 = (b/a) \Jzsoc - x 2 (25) 

dy/dx = b 2 (a - xj/a 2 ? (26) 

Substituting the values at point P in Eq. (25), we have: 



= 28 = (b/a) ^2a(40) - (40) 2 



or: 



b/a = 28/ \/2a(40) - (40) 2 (97) 

and substituting in Eq. (26): 

dy/dx - 0.3 . .ft 2 )^ - 40 > 
(a S )28 



or: 



b/a= ^50 (98) 

fa - 40)/(28) 
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(a): 



Tbese two equations for (b/a) are no. equated in order to solve for 

( *8)/ Jm^Tw? » /c^w ^^^v^ 

(28) 2 ((a - 40)/28) - (0.3) [^(40) - (40) 2 ] 
a = 160 
mtlng this vaiue or (a) oac* i. *. (97), «*.'<«), «- «*» ° f . 

(b) is found to be: 

/ b - (160) (28)// (2) (160) (40) - (40)*/ 

. 42.332 
Tbe equation for tbe ellipse fro, point to point P is then, fro. 

Eq- ( 25 > ! ;.'- ,_ T ..(99)' 

y - (42.552/160) \,520(x) - x* 

w „ have now three equations: the eilipse *- station to Station ' 
«. the supp^entar, — *- station 40 to station ,00, an, ftoa there on 
*. ori g i*al curve. The curvatures T each curve shouid he ne* *<^ T_ 
plo tted, aud a chec* »ade Aether or not the, ,atch ciose*, aad are c OT ti ™ 
XT not, the siope or odinate at station 40 should he chen g ed, and ne. options 
aeveloped until the desired carve continuity is ohtained. 

4 £0 - JU2! * iM t i < fflJ! U* !E ^^^ i • 

* p o«r curve .nations to deveiop oar.ee ,or cross sections - a ho* entai - 
». S a*e procednres as have heen outiined previous*. However, to f a»i^ri 
tt e reader with the shape or the carves „nich can he ohtained, and to rnrthe 
de „onstrate the use o f the t„o tern Power carve action chart or ««. «. 
f e„ motive espies ^ he g iven. T~ ter» actions as derived fto 
ng . „ „»1 he elusive* used, aitho.h it is entire* possiUe, * ^ 
«... to deveiop and use a sin g Xe t- option, or a two or -ore ten* „ 
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developed by the trial and error method of Section 2.60. 

Let us assume that, in Fig. SO, it is required to develop a portion of 
a cross section so that the line starts at point and passes through point P, 
at which (x) = 20, (y) - 20, and the slope is zero. 




FIGURE 50 



Now, obviously, a single term equation could not be developed in this 
case, since it is impossible for a single term equation to have a zero slope at 
any point. Hence, the chart of Fig. 17 can readily be used, as it has been 
developed as a function of the maximum ordinate, at which point the slope is 
zero. 

Suppose, for a first try, that the first and second term exponents are 
assumed to be 0.5 and 0.8, respectively. The equation so far is then: 



y = kjfr) 



°- 5 - yx) ' 8 



.(100) 



From Fig. 17, the maximum ordinate would be at 20.7^ of the chord. For 
the purpose of determining the fundamental curve shape, and solving for (k.) 
and (kg), two points must be evaluated. If the chord is assumed to be 100 units, 
at which (y) = 0, then Eq. (100 ) can be evaluated for the maximum ordinate 
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position and the end of the chord. The «!». ordinate of 20 then occurs at 
(x) - 20.7. Hence: 

y 

and: 



= 20 - tl ( 2 0.7) - 5 - k 2 (20.7)°- 8 - 0»> 



7 



y =: o ^(ioo) - 5 - k,(ioo) - 8 ( 102 ) 

Solving Eqs. (101) and (102) simultaneously: 

kg * 2.95224 
=: N(0.46720) 

and: 

k x = 11.67554 

- N (1.06720) 

Eq. (100) then becomes: 

„ N(l.06720)00 - 5 - N(0.46720)(x) - 8 .........(105) 

How, if the first term exponent is assumed to be 0.5 and that of the 
second, term to be 1.0, going through the same procedure results in the equation: 

y » N(0.90509)(x) ' 5 - tfO-WSOT-DM 1 " — " '"^^ 

. ^tfa the maximum ordinate occuring at 25% of the chord. For first and 
second term exponents of 0.5 and 2.0, respectively, the maximum ordinate occurs 
at 59.6% of the chord, and the equation is: 

y = N(0.62666)(x) - 5 - N(0.62666-5)(x) 2 -° ■••• < 105 > 

With first and second term exponents of 0.5 and 5.0, respectively, the 
ordinate is at 49% of the chord, and the equation is: 

y - N(0.55580)(x)°' 5 - N(0.55580-5)(x) 3 -° ....-(106) 



maximum 
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Changing the exponent of the first term to 0.4, and that of the 
second term to 1.0, the maximum ordinate is at 21.7* of the chord, and the 
equation becomes: 

7 = N(0. 98817) (x) ' 4 - N(0. 78817-1) (x) 1 * (107) 

Finally, increasing the exponent of the first term to 0.6, and leaving 
the second equal to 1.0, the maximum ordinate occurs at 2.8% of the chord, and 
the equation is then: 

y = N(0.83174) (x) * 6 - N(0.03173) (x) 1 * (108) 

Eqs. (103) through (108) are plotted on Fig. 51. To get a more 
comprehensive comparison, however, it is necessary to change the chord in such a 
way that the maximum ordinate position in each case falls on station (x) = 20. 
Each one of the chords will have to be foreshortened in proportion to its 
maximum ordinate position, the ordinates remaining the same while the abscissae 
are changed in the required proportion. Thus each value of (x) in each equation 
is multiplied by a factor equal to twenty divided by the value of (x) at the 
maximum ordinate position. For example, in Eq. (105), this factor will be equal 
to 20/20.7, or 0.9662, and in Eq. (104) the factor is 20/25, or 0.80. 

This has been done in Fig. 32, where all curves are plotted together 
to give a direct comparison. It can clearly be seen that, in the equations 
with an exponent of 0.5 in the first term, the one with the lowest exponent in 
the second term gives the more bulgy curve, while the sections nearest the 
origin become more slender as the power of the second term is increased. It will 
also be noticed, by comparison of Eqs. (104), (107), and (108), that if the 
exponent of the first term is decreased, the curve bulges more near the origin 
even if the exponent of the second term is maintained constant. 
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^ 



Eq. (105) 




FIGURE 51 
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FIGURE 32 



The above method outlines the simple manner in which two term power 
curve equations derived from Fig. 17 can be used for cross section shapes. As 
mentioned previously, single term equations can also be used, provided a zero 
slope is not desired, or two or more term equations can be developed by the 
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trial and error method of Section 2.60. Fig. 32 will serve as a guide, however, 
to the type of shapes obtainable through the use of two term equations derived 
from Fig. 17. 

As a further illustrative example in the use of Fig. 17, suppose that 
a curve is to be developed to pass through points P and Q in Fig. 33, the slope 
at point Q being zero. 
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FIGURE 33 



For a first assumption, let the first term exponent be 0.3, and the 
second term exponent be unity. Then, from Fig. 17, the maximum ordinate of a 
two term equation would be at 17.9/5 of the chord. Since, In this case, this 
maximum ordinate position must coincide with (x) - 20, the chord is: 

- (20)(100)/(17.90) 

- 111.69 

Then, when (y) = 30 at (x) = 20, and (y) = at (x) = 111.69, the 
equations from which (k-^) and (kg) can be found are: 



98. 

k^O) * 3 - k g (20) = 30 

and: 

k-^lll. 69) 0,3 - kg(lll.69) =0 

Solving these two equations simultaneously, we get: 

k 2 » 0.64285 

= N(0. 80811-1) 
and: 

^ = 17.4467 

= N(l. 24163) 

The equation of the curve is then: 

y = N(l. 24163) (x)°* 3 - N(0. 80811-1) (x) (109) 

Evaluating Eq. (109) for (x) = 5, we find that (y) - 25.06 instead of 
the required value of 26.5. If, for a second assumption, the first term exponent 
is left at 0.3, and the second terra exponent is reduced to 0.6, the equation 
becomes: 

y = N(l. 38798) (x) ' 5 - N(0. 69644) (x) ' 6 (110) 

■- For (x) = 5, the value of (y) in Eq. (110) is found to be 26.547, 
which will be considered a close enough check. 

Here again, the trial and error method of Section 2.60 could have been 
used, but very often it will be found that considerable time can be saved by 
the use of Fig. 17, as was done in the illustrative example. The use of either 
method is optional, however, and one will soon find that practice will dictate 
the method most often used. 

4.30 - Transition Curves for Convergent Lines :- For developing a curved 
transition for convergent lines, circular arcs are sometimes employed. While 
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these are easily computed and constructed, they are very often undesirable 
„here it is desired to maintain a smooth fluid flow over the surface, or in 
highway or railroad curves. This is because the straight lines to be joined 
obviously have zero curvature, while circular arcs have constant curvature, 
and therefore the transition from one to the other becomes abrupt. 

To avoid this, power curve equations can be developed that start with 
ze ro curvature, gradually increase to a maximum value of curvature, and then 
decrease in the same manner, finally ending up with zero curvature. Multi-term 
equations can be developed for this type of application by the method given in 
Section 2.60, but the process is laborious, and, in this case, is not needed. 
It is easier to develop a single term equation to fit the curve up to its mid- 
point, and then plot the curve in reverse for the other half of its length. 
This is particularly desirable for cases where there is reverse flow over, or 
on, the surface, such as highways, railroad curves, etc. 

To develop such a single term power curve equation, the problem consists 
basically of finding the equation of a curve that has zero curvature at its 
origin. Let Fig. 54 represent the problem, with the dotted curve being the 
line for which the equation is desirad, and M and NO the lines to be joined. 




FIGURE 34 
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On th8 bisector of the angle 9, point R represents the mid-point of 
the curve. Then the equation for the first half of the curve starts at point 
M and ends at point R. The slope at point R, where the required curve is 
parallel to the line MO, is fixed, and is: 

t = (b)/(a + c) (Ill) 

However, the lengths A and B are not known, and therefore an infinite 
number of curves could be drawn between points M and R. To fix one of these 
unknowns, let us first consider the expression for curvature as given in Eq. (93): 

X - (dVdx 2 )/[l + (dy/dx) 2 ] 1 ' 5 ..(95) 

Since the curvature is required to be zero at the beginning of the curve, 
it is apparent that it will only be necessary to find an equation whose second 
derivative starts with zero and then gradually increases. The value of the ex- 
ponent, as has been demonstrated previously, determines the amount of curvature. 
However, to begin with, it is obvious that an exponent of unity will not fulfil 
the requirements, as this exponent defines a straight line. For an exponent of 
2, there results: 

y = k(x) 2 

dy/dx » 2kx 

d2y/dx 2 = 2k 

Here the second derivative is a constant, meaning a constant curvature, 
which ia also unsatisfactory. Furthermore, any exponent between 1 and 2 will 
have a negative second derivative, indicating a curvature that decreases with 
(x) rather than increasing with it. Hence, the exponent must be larger than 2, 
and then the curve will start with zero curvature and increase in curvature a3 
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W increases, me -**-. °* «- -P-* «"— "*» *" dlStmCe W ' t 
te „«. 34, and the an g le of convene,, ana cen host he «— * * * - «* 

computations. 

With the value of the exponent, (n), assumed, the values of A and B 
« he easi* found. *. *. (=4) in Section *.», for the case of a curve 

v- a * .1th abscissa (A), ordinate (B), and slope (t) at the 
passing through a point with abscissa uw, 

point, we have: :.. ^ 

-. . „ ,/ n = At/B **""*"** 



or: 



B = At/n 



(112) 



* Fig. 54, the triangles ESN and MPO are similar, and hence the slope 
(t) at point R is also given by: 

t = (a - A)/B 



or: 



B = (a - A)/t 



,(H5) 



Equating Eqs. (IK) and (115): 

At/n» (a- A)/t 



or: 



A= (a)(n)/(t 2 + n) 



.(114) 



Haac, the value of A can he found fro. *. (114), and the value of B 
ft. *. (11,). Then the value of (*) for the needed action is siven h 7 *. 
(«) men the option has heea developed, the curvature should he checked to 
^sure that it has not reached a »asi»u» hefore reaching point K, etherise 

* " ' ^ -rpvprsals will occur. Should the maximum _. 
discontinuities in curvature, or reversals, wxxx 

+ ™vmt R a new'value of the exponent has to be assumed, 
curvature not occur at point R, a new v ^ e 
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When plotting the entire curve from points M to 0, it must be remembered 
that, for the first half of the curve, the line HN represents the X axis, with 
point M as the point of origin. For the second half of the curve, the line NO 
represents the X axis, with point as the origin. 

For example, consider Fig. 35, where it is desired to develop the equat- 
ion of a curve connecting points M and 0, with zero curvature at these points 
with respect to the lines MN and NO. Other specified dimensions are: 

a = 200 

b - 132.288 

c = 150 




From Eq. (Ill): 



FIGURE 55 



t * b/(a 4- c) 
- 132.288/550 
= 0.37796 
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Assuming a value of (n) = 3.0, the value of (A) is given by Eq. (114): 
A= (a)(n)/(t 2 +n) ' ■ 

= (200)(3)/[(0.57796) 2 + 3j 



From Eq. (112): 



From Eq. (55): 



190.91 

B ■- At/n 

= (190.91) (0.37796)/(3) 
=24.052 

k = B/(A) n 



- 24.052/(l90.91) 5 
= N [(log 24.052) - 3(log 190.91)] 
.'''■' =N(0. 53867-6) 
The equation of the curve then becomes: 

_ y = k(x) n 

= N(0.53867-6)(x) 3 ...(115) 

The ordiriates, derivatives, and curvatures of Eq. (115) are tabulated 
in Table XIV. . 

It can be seen that this curve will be satisfactory in that the curvat- 
ure continues to increase with (x), and is largest at (x) = 190.91, correspond- 
ing to point R. The minimum radius of curvature, which occurs at point R, is: 
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TABLE XIV 



X 


7 


dy/dx 


dVdx 2 


K 


50.00 


0.452 


0.02593 


0.0010571 


0.001056 


90.00 


2.520 


0.08401 


0.0018668 


0.001849 


150.00 


7.595 


0.17527 


0.0026965 


0.002578 


160.00 


14.160 


0.26550 


0.0033187 


0.002995 


180.00 


20.161 


0.33602 


0.0037336 


0.003181 


190.00 


23.712 


0.37439 


0.0039410 


0.003236 


190.91 


24.053 


0.37798 


0.0041484 


0.003394 



R - l/K 

= 1/0.003394 
- 294.637 

If an exponent of 2.5 had been assumed to fit the case of lig. 55, (a) 
becomes 189.19, (B) becomes 28.603, and the equation Is: 

y = N(0.76418-5)(x) 2 - 5 ..(116) 

The ordinates, derivatives, and curvatures of Eq. (116) are given in 
Table XV below: 

TABLE XV 



X 


7 


dy/dx 


d^/dx 2 


K 


50.00 


1.027 


0.05135 


0.001541 


0.001536 


90.00 


4.465 


0.12401 


0.002067 


0.002020 


130.00 


11.177 


0.21530 


0.002484 


0.002321 


160.00 


18.814 


0.29397 


0.002756 


0.002434 


180.00 


25.256 


0.35077 


0.002923 


0.002456 


189.19 


28.603 


0.37797 


0.002997 


0.002453 
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In this case, it will be noticed that the maximum curvature does not 
quite occur at point R, (x = 189.19), but is instead at (x) = 180. The minimum 
radius of the curve is: 

R = 1/0.002456 
•= 407.6 

For an assumed exponent of 4.0, there results: 

As 193.10 
B = 18.247 

y = N(0.11802-8)(x) 4 (117) 

Table XVI shows the ordinates, derivatives, and curvatures of Eq. (U7): 
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X 


7 


dy/dx 


dVdx 2 


K 


50.00 


0.082 


0.00656 


0.000394 


0.000394 


90. CO 


0.861 


0.03826 


0.001275 


0.001272 


150.00 


3.748 


0.11532 


0.002661 


0.002609 


160.00 


8.600 


0.21500 


0.004031 


0.003767 


180.00 


13.775 


0.30612 


0.005102 


0.004461 


190.00 


17.101 


0.36003 


0.005685 


0.004735 


195.10 


18.247 


0.57796 


0.005872 


0.004806 



The maximum curvature now again occurs at point R, as is desired, and 
the radius of curvature at this point, which is the minimum radius along the 

curve, is: 

R - 1/0.004806 

■ *. 

m 208.07 
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The curves defined by Eqs. (115), (116), and (117) are plotted on 
Jig. (36), and the values of curvature for the three equations are shown on 
Fig. (37): 



7i 




160 



180 



FIGURE 56 




130 



160 180 



FIGURE 37 



By inspecting the curves of Fig. 37, it becomes at once apparent that 
the curve of Eq. (115), for which the exponent is 3.0, is by far the best, since 
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it is very close to a straight line, meaning a uniform rate of change of 
curvature. The equation with an exponent of 4.0, (Eq. 217), is preferred to 
that with an exponent of 2.5, (Eq. 116), since the increase in curvature is 
much less at the beginning, and also because the maximum curvature of Eq. (116) 
actually occurs a short distance before the mid-point of the curve. However, 
the minimum radius with an exponent of 4.0 is smaller than with an exponent of 
2.5, and this may be a point for consideration in some cases. 

The above discussion has been based on the procedure of developing the 
curve so that it is exactly symmetrical about its mid-point, and using different 
reference lines for each segment of the curve. As mentioned previously, such a 
development can be used where the flow over, or on, the surface may be in 
either direction, as in highways, etc. If the flow is in only one direction, 
however, as is usually the case when dealing with the flow of fluids, the 
symmetry of cnirvature about the mid-point may be of lesser importance, and the 
curve may be developed in another manner which is sometimes easier. 

In this method, » and (c) are made equal, and drawing a perpendicular 
through point N will fix (A) equal to (a), as shown in Fig. 38. 




FIGURE 58 
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The slope at point R is again: 

t = b/(a + c) (HI) 

With the value of (A) known, the only unknown is (B), which can be 
found from Eq. (112), and (k) is found in the usual manner from Eq. (53). After 
the equation has been' developed, and the ordinates computed, the curve up to 
point R is plotted with line MP as the reference line and point M as the origin. 

For the second part of the curve between points R and 0, which is plott- 
ed with the same reference line but with point P as the origin, the equation 
developed for the curve MR can also be used providing corrections are made for 
adding in the ordinates of the line NO. These corrections are needed since the 
ordinates which would be obtained from the same equation as for line MR would 
be measured from line NO, and it is desired, for the sake of ease of plotting, 
to use line NP as the reference line. The value of the ordinates of line NO 
is a function of (x), as measured from point P, and is found as follows: . 

\ tan(«0 = b/c 

or: 

b = (c)tan(ot) 

At any point (x) : 

y = (c - x)tan(ot) 

= (c - x)(b/c) 
.- b - (b/c)(x) ....(IIS) 

The required equation for the curve R0 is then: 

y - k(x) n + [b - (b/c)(x)] (119)' 

The exponent (n) is assumed again, as in the first case where the 
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actual mid-point of the curve is used. After the equation has been developed, 
the curvature should be checked to insure that the maximum value occurs at 
point R. If such should not be the case, a new value for the exponent should 
be assumed, and the development repeated. 

As an illustrative example of this procedure, consider Fig. 38 and 
assume the following to be given: 

a - c -A = 100 
b * 44 

The slope at point R is then, from Eq. (Ill): 

t = b/(a + c) 
= 44/(100 + 100) 
= 0.22 

Assuming an exponent of 3.0, and using Eq. (112): 

B = At/n 

= (100)(0.22)/(3) 
= 7.333 
and, from Eq. (53): 

k=B/(A) n 



! = 7.333/(l00) 3 



= N(0.86530-6) 
The equation for the curve MR is now: 

y - N(0. 86530-6) (x) 3 . .....(120) 

The first and second derivatives of Eq. (120) are: 

dy/dx = N(0.34242-5)(x) 2 (121) 
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and: 



d^/dx 2 = N(0. 64345-5) (x) 



(122) 



The ordinates, derivatives, and the curvature of Eq. (120) are given 
in Table XVII below: 

TABLE XVII 



30 



60 



80 



90 



0.0 



0.198 



1.584 



3.755 



dy/dx 



0.0 



0.0198 



0.0792 



0.1408 



5.346 



100 



7.333 



0.1782 



dV«3bc S 



0.0 



0.00132 



0.00264 



0.00352 



0.2200 



0.00396 



0.0 



0.00132 



0.00262 



0.00542 



0.00578 



0.00440 



0.00410 



For the second portion of the curve, R0, we get, from Eq. (119): 

" ■? 

y = N(0.86530-6)(x) 3 * [44 - 0.44(x)] .....(123) 



and the first and second derivatives are: 
dy/dx = N(0.54242-5)(x) - 0.44 



.(124) 



and: 



d 2y/<3x 2 = N(0.64345-5)(x) 



.(125) 



Table XVIII lists the ordinates, derivatives, and curvatures as obta^ed 
from Eqs. (123), (124), and (125), for the curve BO. 

While the curvature for the curve RO starts with zero, and ends up at 
point R with a value of 0.00410, the same as for the curve MR, the intermediate 
values are slightly greater due to the somewhat shorter length of the second 
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TABLE XVIII 



X 


7 


dy/dx 


dV/dx 2 


K 





44.000 


-0.4400 


0.0 


0.0 


30 


30.998 


-0.4202 


0.00132 


0.00103 


60 


19.184 


-0.3608 


0.00264 


0.00220 


80 


12.555 


-0.2992 


0.00352 


0.00309 


90 


9.746 


-0.2618 


0.00396 


0.00559 


100 


7.333 


-0.2200 


0.00440 


0.00410 



.('. 



part of the curve. 

The values of curvature should now be plotted to determine their 
general shape. After doing so, it may be found that the rate of change of 
curvature may be improved by the assumption of another value for the exponent. 
It will be found in all cases that it is well worth the time to plot the 
curvatures, as sometimes undesirable characteristics will be brought out that 
are not readily apparent from tabular data . 

4.40 - Transition Curves for Parallel Lines ;- In the previous section 
the development of transition curves for convergent lines was demonstrated. 
For parallel lines, the procedure is the same after a diagonal has been drawn 
through the parallel lines in such a way that two sets of identical convergent 
lines are obtained. However, before an equation can be developed, it is first 
necessary to find a line, or lines, that will divide the system into two ident- 
ical patterns without an excessive slope. To determine this line, let point M 
and point T be the points of tangency, as shown in Fig. 39, and let (a) be the 
length, and (b) the distance between, the two parallel lines. 

Then the line PE can be drawn perpendicular to the two lines to be join- 
ed, MU and TIT, so that MP - PU - (a/2). Point S can then be so located on line 
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PE that SP = SE = (b/2). The problem is now to locate line NSV so that all 
reference lines, MN, NS, SV, and VT, are equal. To do this, the diagonal line 
MST is first drawn connecting points M and T, and forming an angle ^ with 
lines BUT and FT. The length of the line MS is: 



MS - ]/(W) Z + (PS) 2 

and the angle ^ is: 

tan(^) = SP/MP 

Considering now the triangle MFN, where M is a perpendicular to line 
MS drawn through the mid-point of line MS, it is seen that: 

cos(/) - MF/MN 

and therefore the required length MN is: 

MN = (MF)cos(^) C 126 ) 



Lines MN, NS, SV, and VT will now be of equal length, and, if we now 
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consider the convergent lines MNS and SVT, it is readily seen that each of 
these is the same case as that of the convergent lines considered in Fig. 35, 
Section 4.50. Hence, the equation for the curve MG may be found by the first 
method outlined in Section 4.30, (when using the different reference lines of 
equal lengths), the slope at point G being equal to (b/a). This curve is used 
four times, as each quarter will use the same curve. The first point of origin 
will be at point M, with reference line MN, the second point of origin will be 
at point S, with NS as the reference line, and so forth. 

The problem may also be solved by use of the second method given in 
Section 4.30, where point G lies on the perpendicular to line MP drawn through 
the mid-point of line MP, and only one reference line is used. Once the equat- 
ion for the first section of the curve, MG, has been determined in the form 
y - k(x) n , the equation for the second quarter becomes: 

y = k(x) n + [(Tj/2) - (2b/a)(x)] (127) 

For the third quarter, the equation is: 

■ * 

y - -k(x) n + [(b/2) + (Zb/a)(r)] (128) 

and for the fourth quarter: 

y = b - k(x) n (129) 

All curves are plotted from the reference line HIT, and have respective 
origins at points M, S, S, and T. 

4.50 - Transition Curves fo r Divergent Lines:- Transition curves for 
divergent lines can also be developed in a similar manner as demonstrated in 
Section 4.30 for convergent lines, either by using different reference lines 
of equal length, as shown in Fig. 34, or by the n mid-span n method, with one 
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reference line, of Fig. 38. 

Consider first the method of different reference lines, as shown in 
Fig. 40, where lines MS and QG are to be joined by a smooth curve tangent at 
points M and G. It can be seen that the whole problem actually resolves itself 




FIGURE 40 



,1 



into finding a diagonal line so that the distances MN, NO, 0Q, and QG are 
equal, since, once this line has been found, the equation for the curve for 
each quarter can be found exactly as outlined in Section 4.50. It should be 
noted, however, that the equations for the curve from points M to 0, and from 
points to G, will Sot be the same in all cases, as point will not 
lie on a straight line through points M and G. 
If the distance MN be called (a), then: 

MN a NO = 0Q = QG = (a) 

The following equalities can be set up by reference to Fig. 40: 

QV = (QG)sin(e) 
= (a)sin(e) 
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NT = /(N0+0Q) 2 - (QT) 2 

- |/(2a) 2 - (b + (a)sin(e)) 2 
= V4a - (b + (a)sin(e)) 

VR - YS - (a)cos(e) 
Therefore: 

d » (a) + (a)cos(e) + 74a 2 - (b + (a)sln(e)) 2 (130) 

Bq. (150) may be re-arranged as follows: 
[(1 + cos(e)) 2 + sin 2 (e) - 4> 2 + [(2b)sin(e) - 2d(l + cos(e))]x + b + d '- 0..(131)> 

This is a quadratic equation which may be solved by means of the standard, 
quadratic formula: 

a = -E±l/E 2 -4DF/(2D) ........(132) V 



where: 

D 



[(1 >cos(e)) 2 + sin 2 (e) -4] ....... ...:(l32a) 

E »[(2b)sin(e) - 2d(l + cos(e))] .,..(132b) 

F = (b 2 + d 2 ) (132c) 

Also: 

sin(e) = QT/NQ - : --" 

= (b + (a)sin(e))/2a -(133) 

h- (a)sin(e) -(134) 

Thus the value of (a) needed to form four equal reference lines, MN, 
NO, 0Q, and QG, can be found from Eq. (132), and the values of (6) and (h) 
found from Eqs. (133) and (134), respectively. Then, by drawing the line NR 



M 
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to bisect the angle MNO, the mid-point of the curve is located at point R. 
The slope of the curve at point R is: 

t = h/(a + c) 

and the values of (a), (B), and (k) can be found, as outlined in 
Section 4.30, by means of Eqs. (114), (112), and (53), respectively, after 
assuming the value of the exponent (n). The equation can then be put in the 
usual 7 = k(x) n form. 

As an example in the use of this method, assume that the following are 
given with respect to Fig. 40: 

d » 850 
b = 200 
e - 10* 

It is then first necessary to find the required value of (a) as given 
by Eq. (132). Hence, from Eqs. (132a), (132b), and (132c), we have: 

D » [(1 + cos(e)) 2 + sin 2 (e) - 4] 
= [(1.9848) 2 + (0.1756) 2 - 4] 
= -0.0304 

E =[(2b)sin(e) - (2d)(l + cos(e))] 
= [(400) (0.1736) - (1700) (1.9848)1 

- -5504.72 

F = (b 2 + d 2 ) 

- (200) 2 + (850) 2 
= 762,500 
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Then: 



a = 



-(-5504.72) ± \/(-5504.72) 2 - 4(-0.0504) (762500) 

(2) (-0.0504) 



= 230.243 



From Eq. (133): 



sxn 



(0) - (200) + (250.245) (0.1756) 
~ : (2) (230.243) 



- 0.52115 
From, trigonometric tables, corresponding to the above value of sin (9): 

cos(8) = 0.85547 
tan(6) = 0.61062 

Also, referring to Fig. 40: 

m - (a) cos (9) 

= (250.243) (0.85347) 

= 196.505 
m = 196.505 

IS = (a)cos(e) 

= (230.243) (0.98481) 
= 226.746 

h '■ (a)sin(e) 

- (250.245) (0.52115) 
= 119.991 

These dimensions are shown on Fig. 41 for reference. 

We may now proceed with the development of the equation for the curve 
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between points H and 0. The slope at point R, (Fig. 40), is given by Eq. (Ill), 
except that, in this case, (b) in Eq. (Ill) is represented by (h). Hence: 

t ■ h/(a + c) 
;- (H9.99l)/(426.748) 

=0.28118 

. ■ ■ ''"'■■■■' 

Assuming an exponent of 3.0, the value of (A) is found from Eq. (114) r 
A= »(n)/[(t)% (n)] 



(230.243) (3)/ [(0.28118) 2 + 3] 



s 224.331 



From Eq. (112) ! 



and, from Eq. (53): 



B = At/n 

- (224.531) (0.28118)/(3) 
a 21.026 

k = B/(A) n 



-A 



it) 
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= N[(log 21.026) - 3 (log 224.331)] 
» N(7. 32276-6 - 7.05267) 
= N(0.27009-6) 

The equation for the curve between points M and R, and points R and 

0, is then: 

y = N(0.27009-6)(x) 5 (155) 

This curve is plotted with the reference line MN, and the point of 
origin at M. For the curve between points and G, Fig. 40, it is first necessary 
to find the slope of the line 0G with respect to the new reference line 0Q. The 
slope of the line 0Q has already been found: 

sin(e) =0.52115 
(9) - 30°24 t 34 n 



Also: 



tan(y) = (b - h)/(W + YS) 

= (200.00 - 119.991)/(196.505 + 226.746) 
=0.18903 
(r) =10°53 , 47 n 



Hence the angle between lines OG and 0Q is: 



= (e) - (r) 



20*30 f 57 n 



and the slope of line 0G with respect to the reference line 0Q is: 

t = tan(20*30 , 57 B ) 
= 0.37420 
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Proceeding as in the case for the curve between points M and 0, and 
employing Eqs. (114), (112), and (53) r 

A - (250.243) (3)/ [(0.37420) 2 + (3)] 

* 219.975 
B = (219.975) (0.37420)/(3) 

= 27.458 
k = N[log(27.458) - 3 log(219.975)] 

= N(0.41124-6) 

and the equation of the second curve is: 

y - N(0.41124-6)(x) 3 ..(156) 

It must be remembered that Eq. (135) is used twice to complete the 
curve from point M to point 0, the first segment being between points M and R, 
with the origin at point M, and the second segment between points R and 0, with 
the origin at point 0. Similarly, Eq. (136) is used twice to plot the curve 0G, 
point being the origin of the first segment, and point G the origin of the 
second segment. 

Tables XIX" and XX show the calculations from Eqs. (155) and (136), 
respectively. The values are given for only one segment of each curve, as quite 
obviously the second segments in each case are the same, with respect to their 
reference lines, as the first segments. 

The transition curve for divergent lines can also be plotted from one 
reference lin8, instead of the multiple reference line method outlined above. 
To do this, however, the equations must be developed in a different manner. 
Consider lig. 42, for example, and let it be supposed that the transition curve 
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TABLE XIX 



X 


7 


dy/dx 


dVdx 2 


K 


0.0 


0.0 


0.0 


0.0 


0.0 


100.00 


1.862 


0.05587 


0.001117 


0.00111 


180.00 


10.861 


0.18105 


0.002011 


0.00192 


220.00 


19.830 


0.27042 


0.002451 


0.00221 


224. 33 


21.025 


0.28118 


0.002507 


0.00224 







TABLE XX 






X 


7 


dy/dx 


d^/dx 2 


K 


0.0 


0.0 


0.0 


0.0 


0.0 


100.00 


2.578 


0.07733 


0.001546 


0.00153 


180.00 


15.033 


0.25055 


0.002784 


0.00254 


210.00 


23.873 


0.34104 


0.003248 


0.00275 


219.975 


27.438 


0.37420 


0.003402 


0.00280 




FIGURE 42 
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indicated by the dotted line, MG, is to be developed so that it can all be 
plotted from the one reference line MS. 

To do this, the distance (d) is divided into four equal parts by the 
line NQ. The slope of this line is then: 

tan(6) = [tan(e)(d/4) +b]/(d/2) (137) 

The distance (h), which locates the intersection of the diagonal line 
NQ with a perpendicular at distance (d/2) from the origin of the curve at point 

M, is: 

h = tan(e)(dA) .(138) 

The. slope of the curve UK) at its mid-point i3 then: 

t =h/(d/2) 

Hence, having determined these quantities, and assumed an exponent, the 
equation for the first quarter of the curve, MR, can be found in the usual 

manner from:-. 

A= (d/4) 

B = At/n .....(112) 

k = B/(A) n ..........(53) 

and: 

y = k(x) n , r . 

For the second quarter of the curve, RO, it is only necessary to add 
to the standard curve equation the ordinates of the line NO, since, it will be 
remembered, the ordinates of this section of the curve with respect to the line 
NO are the same as the ordinates of the first quarter of the curve with respect 
to line MN. The origin, in this case, is at point W, the distances (x) being 
measured to the left from this point. Hence, the curve equation is: 



■-M 



W 
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y = [h - tan(6)x] + k(x) n (139) 

For the third portion of the curve, OE, a new equation must be set up, 
since the slope of the curve at its mid-point, E, is different from that of the 
curve MO. This slope is equal to the slope of the line OG, and is, with respect 

to line MS: 

...V = (b - h)/(d/2) ..........(159a) 

Since the value of the slope of the curve at its mid-point used in 
Eq. (112) for finding (B) is given with respect to the basic reference line of 
the curve, (which is line OQ in this case), the value of (^) represents this 

slope, and is: 

tanOzO = tan(6) - t' ...............(139b) 

Since (A) is still equal to (d/4), the value of (B*) with respect to 

line OQ is: 

B* - A tanOzO/n .. ...... ............... ........(140) 

The ordinates of line OQ, measured from the reference line MS, are 

given by: 

y = h + tan(8)x 

where the origin is at point W and the distances (x) are measured to 
the right from this origin. Using the value of (B*) from Eq. (140) in Eq. (53) 
to find (k'), the equation for the third portion of the curve, with respect to 

* 

the reference line MS, becomes: 

y = h + tan(0)x-k'(x) n .(141) 



• 



The 



origin in this case is at point W. For the curve EG, the equation 
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will be: 

y = b + tan(e)x - k*(x) n ( 142 > 

where the origin is now at point S, and the distances (x) are measured 
to the left from this point. Here again, in finding the value of (B) in Eq. 
(112) to use in determining (k' ) for Eq. (142), care must be taken to refer the 
slope of line OG to the reference line QG about which the basic equation of the 
curve is developed. The slope, (O, is therefore: 

t n =t* + tan(e) - — " < 142a ) 

and therefore: 

t ■«/ (142b) 

B = At /n • v 

For an example of this procedure, assume the following are given in 

Fig. 42: ' 

d = 850 

b = 200 

e = 10° (tan(e) = 0.17633) 

n = 3.0 (assumed) 

The slope of the diagonal dividing line, SQ, is, from Eq. (137): 
tan(6) =[(0.17633)(212.5) + (200)]/(425) 

- 0.55875 
and, from Eq. (138): 

h - (0.55875) (212.5) 
= 118.735 

The # slope of line M0, which is also the slope of the curve M0 at its 
mid-point, is: 
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t = (ll8.735)/(425) 
= 0.27938 

Hence, since (A) = (d/4) = 212.5, (B) is, from Eq. (112): 

B = (212.5) (0.27938 )/(3) 
= 19.789 

From Eq. (53): 

k - (19.789)/(212.5) 5 

. N(0. 31435-6) 
The equation of curve MR is therefore: 

y = N(0.31435-6)(x) 5 .....(143) 

For the curve R0, with its origin at point W, it is only necessary to 
correct the above equation by the amount of the ordinates of line NO. This is 
done in Eq. (139), and- thus the required equation is: 

i 
4 

y = [ll8. 735 - 0.55875(x)] + N(0. 31435-6) (x) 5 (144) 

In developing the equation for the third portion of the curve, OE, it 
is first necessary to determine the slope of the line OG, (which is the slope 
of the curve at its mid-point), with respect to line OQ. This is because the 
curve equation is actually first developed with respect to line OQ in the stand- 
ard manner, and then corrections applied to take into account the ordinates of 
this line. Hence, the slope of line OG is, with respect to line MS: 

t' = (b - h)/(d/2) 

= (200 - 118.735)/425 
= 0.19121 



t«= 
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From Eq. (139b), the slope can be corrected to be with reference to 
line OQ: 

tan(frf) = 0.55875 - 0.19121 
- 0.36754 
and therefore: 

B f = (212.5) (0.36754 )/(3) 
= 26.034 ' 



and: 



k' = (26.034)/(212.5) 3 



= N(0.43546-6) 
The equation for the third portion of the curve is then, from Eq. (141): 

y -[ll8.735 + 0.55875(x)] - N(0. 43346-6) (x) 3 (145) 

For the fourth, and last, portion of the curve from point E to point 
G, the same procedure is followed as for the third quarter of the curve. The 
slope of line 0G with respect to line MS is 0.19121, as already found, and 
therefore, referring it to line QG so that the value of (b) may be found, the 
slope becomes, from Eq. (142a): 

t" = 0.19121 + 0.17633 
= 0.36754 

From Eq. (142b): 

b' -(212.5) (0.36754)/(3) 
* 26.034 

and, using Eq. (53): 



U 



\> 
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k' = (26.034)/(212.5) 3 
- N(0.43346-6) 

Eq. (142) then gives the equation for the fourth portion of the curve, 
defined by the line EG: 

y = 200 + 0.17633(x) - N(0.43346-6)(x) 3 (146) 

The equations for plotting the four portions of the curve from the one 
reference line, MS, are then given by Eqs. (143), (144), (145), and (146). The 
ordinates, derivatives, and curvatures for the four portions are tabulated in 

Table XXI. 

From inspection of Table XXI, it will be seen that the ordinates and 
curvatures match at all the junction points, indicating a continuous curve. It 
will be noted, however, that the curvatures for curves MR and HO, and 0E and 
EG, are not exactly the same along the entire length of each respective curve, 
due to the fact that/when using this method of developing the equations, points 
R and E are not quite' at the mid-spans of their respective curves. This differ- 
ence in curvature seldom, if ever, is of any consequence, particularly when 
considering flow on, or over, the developed line in only one direction. However, 
as has been noted before, the curvatures can be made symmetrical about their 
exact mid-points by developing the equations in the manner outlined previously, 
consisting of four separate reference lines. 

4.60 - fl n pli nation of Pow*r 0ur v *« to Projectiles qttd SttWfKmiC Shades :- 
Many bodies designed for high speeds require pointed, straight leading edges 
in order to reduce or delay the effects of compressibility. The lines for these 
bodies can be developed in the same manner as outlined in Section 2.60, except 
that the exponent of the first term becomes unity. The first term therefore, is 
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TABLE XXI 
Ordinates, Derivatives, and Curvatures for Transition Curve of Fig. 42 



Curve WJS. 












X 


7 


dy/dx 


dVd* 2 


K 


0.00 


0.000 


0.00000 


0.000000 


0.00000 


100.00 


2.062 


0.06187 


0.001237 


0.00123 


180.00 


12.027 


0.20045 


0.002227 


0.00210 


200.00 


16.499 


0.24747 


0.002470 


0.00227 


212.50 


19.789 


0.27938 


0.002630 


0.00234 


Curve RO 












0.00 


118.735 


-0.55875 


0.000000 


0.00000 


100.00 


64.922 


-0.49688 


0.001257 


0.00089 


180.00 


50.187 


-0.35830 


0.002227 


0.00186 


200.00 


25.484 


-0.31128 


0.002270 


0.00216 


212.50 


19.789 


-0.27938 


0.002630 


0.00254 


.Curve OE 




> 


0.00 


118.735 


-0.55875 


0.000000 


0.00000 


100.00 


171.897 


-0.47736 


0.001630 


0.00112 


180.00 


203.488 


-0.29504 


0.002930 


0.00259 


200.00 


208.780 


-0.23318 


0.003260 


0.00301 


212.50 


211.436 


-0.19121 


0.003460 


0.00328 


Curve EG 












0.00 


200.000 


-0.17633 


0.000000 


0.00000 


100.00 


214.920 


-0.09494 


0.001650 


0.00143 


180.00 


215.917 


-0.08738 


0.002930 


0.002S2 


200.00 


213.563 


-0.14924 


0.003260 


0.00515 


212.50 


211.436 


-0.19121 


0.003460 


0.00328 
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a slope. 

For instance, considering Fig. 43, assume that the desired leading edge 

angle is 18 degrees. Then, since the exponent is unity, the value of the const- 
ant (k) for the first term is found from Eq. (53), where (b/a) is equal to the 
tangent of the angle, or: 

k = tan(18°) 
* 0.3249 
and, in the usual log form: 

7 - N(0.51175-l)(x) ( 147 > 




FIGURE 45 



Now, if the second point through which the curve must pass is point N, 
with an ordinate of 15, an abscissa of 100, and a slope of zero, the second term 
may be found by the method outlined in Section 2.60 for multi-term equations. 

Substituting in Eq. (147) to find the value of (y) *hen (x) is 100: 

y = 0.3249(100) 
= 32.49 

Since the required ordinate is 15, the correction to be applied to the 
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ordinate of the first term equation is, from Eq. (69): 

C* - E - B* 

= 32.49 - 15.0 
= 17.49 

The first derivative of Eq. (147) is: 

dy/dx - 0.3249 

and, since this is also the slope at (x) = 100 for the first term 
equation, the slope correction to be applied by the second term of the equation 
is, from Eq. (71): 

D' =t-t' 

- 0.3249 - 
= 0.3249 

The required exponent of the second term is found from Eq. (73): 

v z = (aV)/(c») 

= (100) (0.3249 )/(l7.49) 
= 1.858 

and, from Eq. (72): 

k 2 = (cW)* 2 

= (17.49)/(100) 1 * 858 

- N [log(17.49) - 1.858 log(lOO)] 
= N(0. 52679-3) 

Since C and D are positive, the second term is negative, and then the 
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equation becomes: 

y = N(0.51175-l)(x) - $(0.52579-5) (x) 1 ' 858 (148) 

More terms can be added, if desired, in order to develop the curve for 
additional points. This procedure will not be outlined again here, however, as 
this has already been done in Section 2.60. 
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CHAPTER V 

MENSURATION OF THE LINES, AREAS, AND VOLUMES FORMED 
BY POWER CURVE EQUATIONS 

5.00 - Length of a Line Defined by a Power Curve Equation;- In many cases it 
may be desired to determine the length of a line which has been plotted from 
a power curve equation. This can be done mathematical^ in a few cases, but 
even then the equations become unduly complicated. Hence, it is far simpler 
to determine the length graphically, making use of the following equation from 
the Calculus: 

S . /[l ♦ (dy/dr)«]*(dx) ....(149) 

where (S) is the distance from (x) - (a) to (x) - (b), and (dy/dx) 
the slope of the curve at any distance (x) . 

Since we usually have the squares of the derivatives from the curvature 
computations, it is only necessary to add these values to unity, for every stat- 
ion of (x), and plot them against (x). The square root of the area under this 
curve, which can be computed, or obtained by means of a planimeter, will then 
give the length of the line to any desired station. 

5.10 - Area Under a Power Curve ;- The area under any curve developed by a power 
curve equation may be easily determined mathematically by means of integration. 
This may be explained by reference to Fig. 44, which shows a curve plotted from 
the power curve equation: 

y . ^W* - *W"« - ' < 150) 
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jf a very small element is considered, with an infinite ssimal thickness 
(dx), then (y) will be the average height, and the area of the element will be 
(y)(dx). Now, if another element of thickness, (dx), is taken next to the first 
one, and so on until the entire area is taken, the total area will be the sum 
of all the little areas. While (dx) is the same for each element, (y) varies, 
and therefore the total area under the curve will be the sum of all the (y'3) 

multiplied by (dx), or: 

A^/y(dx) 

To determine the area between points F and E, these limits must be 

specified as follows: 

A=A(dx) ----- - ( 150 »> 



Substituting the value of (y) from Eq. (150): 
A^^xfl-Mxf^dx) 

' = /^(xf^dx) -/\(*) n 2(dx) 
T. F 

Integrating the above, we obtain for the area: 



.(150b) 



A = 



kl(x) 



n 1 + l 



IE 



n l + 1 



. .no + 1 
k 2 (x) 2 



n 2 + l 



E 



F 



.(151) 



To solve the above equation, the value of (x) at station E is first 
substituted, and the equation solved for (A), the area from (x) = to (x) - E. 
Then, if the (x) at station F is anything but sero, where (*) is the horizontal 
distance from the origin of the curve, this value of (*) is substituted and 
the equation again solved for (A), where this time the area is from (x)- to 
(x) = F. The final answer, giving the area from (x) - F to (x) - E, is the 
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difference between these two values of (A). Obviously, if the (x) at station 
should actually be zero, no second substitution need be made, and the area 
is that found by substituting only the (x) at station E, since, as noted above, 
this area is the total area from (x) • to (x) = E. 

Eq. (151) applies for a two term equation. For an equation with more 
terms, the procedure is the same as outlined above, except that additional terms 
are added. 

As an illustrative example, let it be supposed that it is desired to 
find the area, from (x) = to (x) = 262 inches, under the curve defined by the 
equation: 

y » N(0.66636)(x) ' 4 - N(0.28565-5)(x) 2 ' 625 

From Eq. (150b), the area is given by: 

y-262 y.262 

A */ N(0.66656)(x)°- 4 (dx) -/ N(0. 28565-5) (x) 2 - 62b (dx) 

integrating the above equation, the area i3 given by Eq. (151): 

t262 



A = 



Hfo.eees sUx) 1,4 



1.4 



N(Q.28565-5Ux) 5 - 625 



5.625 



262 




Substituting (x) = 262 inches in the above equation, we can solve for 
the area from (x) = to (x) - 262 inches as follows: 

log(262) = 2.4183 
log(1.4) = 0.14613 
log(3.625)= 0.55951 

log(A) = [0.66636 + 1.4(2.4183) - 0.14613] - [(0.28565-5) + 3.625(2.4183) - 0.5593l] 
- 3.90585 - 3.49268 



rfi: 



i 
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A = (8051 - 3109) square inches 
= 4942 square inches 



For the lower limit of (x) =0, no substitution need be made, and hence 
the required area is as found above, or 4942 square inches. 

5.20 - Center of Gravity of an Area Under a Power Curve t- In the previous section, 
the area under a curve plotted from a power curve equation was found by summing 
up the areas of all the little elements between certain limits. If now each little 
element is multiplied by its distance from a fixed point, the result is the moment 
of this area about that fixed point. This fixed point, or center of the moment, 
can be assumed to be located anywhere, but for convenience it is preferably taken 
at the origin so that the distance from the center of moment will equal some (x). 

For each element, then, the moment will be (x)(y dx), and, for the entire 
area, the moment will be: 

M-/x(ydx) ............(152) 

But, from Eq. (150a): 

. A -/y(dx) 

Therefore, the moment is equal to the area multiplied by the distance 

(x), or: 

M = A(x) 

or: 

x - H/A ..(153) 

Eq. (153) then gives the expression for the distance from the center 
of moment, or origin, to the center of gravity. 

Assuming, for example, that the equation of the curve is again expressed 
as in Eq. (150), and substituting for (y) in Eq. (152), the moment of the area 
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from (x) m p to (x) « E is: 

M "^pi (x)ni - 1 %W n2 ]wW 



2 + V) 



Integrating the above, we obtain for the moment: 



fcj.(x) 



ni + 2 



TE 



ij + 2 



*2 



/ \*2 + 2 

(x) <■ 



n 2 + 2 



IE 



F 



(154) 



For example, assume that it is desired to find the center of gravity 
of the area under the same curve as illustrated in the previous section, Fig. 44, 
and defined by: 



7 - N(0.66636)(x)°* 4 _ N(0.28565-5)(x) 2 * 625 



Since the area has already been found, it is now only necessary to find 
the moment from Eq. (154): 



M = 



N(0.66656)( x ) 2 * 4 



i262 



2.4 



N(0. 28565-5) (x) 4 - 625 
4.625 



262 



Substituting (x) = 262 in the above, the moment is found as follows: 

log(262) m 2.4183 
log(2.4) =0.38021 
log(4.625) s 0.66511 

g(M) = J0.66636 + 2.4(2.4183) - 0.3802l] - [(0.28565-5) + 4.625(2.4183)-0.665il] 
= 6.09007 - 5.80518 
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M = (1,250,514 - 658,528) 
* 591,986 in 3 

Since th8 area under the curve, from (x) = to (x) = 262, was found 
in the previous section to be 4942 square inches, the distance of the center of 
gravity from (x) = is found from Eq. (153), or: 

x - 591986/4942 
= 119,79 inches 

The above example illustrates the procedure for a two term equation. How- 
ever, the procedure to be followed for three or more term equations is exactly 
the same as outlined above, except for the addition of more terms in the same 
manner as the second term was handled above. 

5.50 - Surface Area of a B ody Tteveloned from a Power Curve;- The surface, or 
wetted, area of a body is sometimes required in the computation of skin friction 
drag, or in the computation of the weight of the body. In the case of a body of 
revolution, it can be determined by simple integration in the following manner: 

Consider the body shown in Fig. 45, which represents a body of revolution 
of which the surface area is to be found. 




FIGURE 45 



138. 

For each small element of surface area, with a width (dx), the area 
ia 2Tf(r)(dx), where (r) is the radius. However, (r) is equal to (y), the ordinate 
of the curve at distance (x) from the origin. Hence, summing up all these element- 
al areas, the total surface area is: 



*■-/'' 



2-rr( 7 )(dx) 

F 



or, since 2Tfis a constant: 



-/< 



W=2y (y)(dx) (155) 

F 

Comparing Eq. (155) with Eq. (150a) in Section 5.10, which gives the area 
under a power curve, it can be seen that the only difference is the multiplier, 
2 IT. The surface area of a body of revolution developed from a power curve equat- 
ion is therefore equal to the area under the power curve, as found by Eq. (150a) 
or (151), multiplied by the constant 2TT. Since the procedure for finding the 
area under the curve has already been outlined and demonstrated in Section 5.10, 
no additional illustrative examples need be given. 

For a body with an elliptical cross-section, where the contour lines in 
plan and elevation views are developed from two different power curve equations, 
the exact calculation of the surface area becomes quite complicated. However, a 
reasonably close approximation can be made, if the difference between the major 
and minor axes is not too great, in the following manner. 

If the plan and elevation view contour lines are developed from two 
equations, designated (y-j_) and (yg), the surface area is approximately: 



lt» 2TT 

F 



yi fa +y*)/ 2 ]< d3C > 



Iff 7l (dx) + Tf/ygCdx) 
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...(156) 



Eq. (156) can then be put in similar form to Eq. (151), the number of 
terms depending on the number of terms in the equations (yj) and (yg). 1* caa 
then be solved as outlined and demonstrated in Section 5.10, remembering, of 
course, to multiply by the constant It . 

5.40 - Volume of a Body of Revolution Developed from a Power Curve t- The body 
of revolution shown in Fig. 46 was formed by rotating the line developed from a 
powBr curve equation 560 degrees about the X axis. Let it now be assumed that 




FIGURE 46 

the internal volume of the body so formed is to be found. 

If a circular disc of infinitessimal thickness, (dx), is considered, 
with a radius equal to (y), the ordinate of the curve at a distance (x) from 
the origin, then the volume of this disc is: 

V = TTy^dx) 

Now, summing up all such discs along the entire length of the body, the 
total volume becomes: 



*"£> 



y*(dx) (157) 



zr 
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If the equation of the power curve equation is of the following form, 
for example, 

7 = k x (x) ni - kg(x) n 2 



then, substituting for (7) in Eq. (157): 
V 



= Trj&xf 1 - * 2 W n2 ] 2 ( fe ) 



.(158) 



The expression in brackets can be squared, as follows: 



r2 '- (k^Cx)**-- ^^(xfl + n 2 + (k/tx)^ .....(158a) 



and, substituting back in Eq. (158): 



(fc L ) 2 (x) 2n l(ebc) - 7y 



* It/ ( 



L 2, x 2n, 



Skik^x) 11 ! " D 2(dx) ^ -jr/ (k2)"(x)^2(dx)...(l59) 



Integrating Eq. (159), we get: 



V = 



2. .2n T + 1 
IT^) (x) J- 



2n 1 + 1 



tE 



jf 



ZTT Ttjk^ix) 1 



n, + n + 1 



nE 



n l + n 2 + 1 



. .2- \2n„ + 1 
TTtkg) (x) 2 



tE 



2n 2 + l 



, . (160) 



jf 



Eq. (160) then gives the expression for the volume of a bod7 of revolut- 
ion formed b7 rotating the curve of a two term power equation 560 degrees about 
the X axis. For basic curve equations with more or less terms than two, Eq. (160) 
will be of the same form except for the addition or subtraction of more terms, 
as determined b7 substitution in the fundamental equation for the volume, Eq. 

(157). 

As an illustrative example, consider the previousl7 demonstrated equation: 



T*~ 



T gtti ^ ir m i r a ¥¥»-*-ft- i ^^ 
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4 ? 625 

y - N(0.66636)(x) ' - N(0. 28565-5) (x) 



Squaring this equation for substitution in Eq. (157), we get, from 
Eq. (158a): 

y2 = N(l.33272)(x) 0,8 - N(0.25304-4)(x) 5 *° 25 + N(0.57130-10)(x) 5 * 25 ° 



If the volume is desired between tb.8 same limits as before, (x) 
(x) = 262, then, from Eq. (159): 



= to 



Ttf N(l.33272)(x)°* 8 (dx) - TXf N(0. 

.262 
+ TlV N(0. 57130-10) (x) 5 * 25 (dx) 



25304-4) (x) 3 - 025 (dx) 



By integrating the above, the volume becomes as given in Eq. (160): 



( TT)N(1. 55272) (x) 



1.8 



1.8 



262 



(TT ) N(0. 25504-4) (x) 



4.025 



1262 



4.025 



(TT )N(0. 57150-10) (x) 



6.25 



6.25 



262 



log(V/TT) = 5.43039 - 9.38193 - 4.88970 
(7/tr ) - 269394 - 240950 - 77587 
= 106031 
V = 333,106 cubic inches 
= 192.77 cubic feet 
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5.50 - Volume of a Body of Elliptical Cross Section Developed from Two Power 
Curve Equations ; - In Section 5.40 the method was demonstrated whereby the 
volume of a fe&dy of revolution could be found where the contour lines were 
defined by a power curve equation. Somewhat the same procedure is followed in 
the case of a body of elliptical cross section, where the contour lines in both 
plan and elevation views are developed from power curve equations. Here, however, 
instead of squaring the one equation for (y), we multiply together the two equat- 
ions in (y)l, and then solve as in Section 5.40. 

For example, assume that the two equations defining the contour lines 
Jin .plan and elevation views are: 

*a = ; V X) ^ " V x) " a2 •• ......(161) 

r b - ■ V 1 '" 1 * - V l)Bhe ' • • " (162) 

Jj The volume is then given by Eq. (163): 

(y a )(y b )dx .(163) 



* m1t / t 



F 
Substituting Eqs. (161) and (162) in the above: 

y^V 1 ^ + nbl(dx) "'y V'i (x,Dbl + n&2(d3c) 



V 

F v F 



-Tr/kb^^x)^ + ^(dx) + 7r/ E k b2 k a2 (x) na 2 + Db 2(dx) (164) 

°F : ' ■ F' 

Eq. (164) then gives the final expression for the volume of the body 
between the arbitrary points F and E. The integration of Eq. (164) is: 






■ : M: 



V » 



H3. 
n__ + n^ + ll E 



t n *l^ ,l bl + 1 * ,B 



**! + n *l * X 



- etc.) etc ..(165) 

F 



As an illustrative example, assume that It is required to find the 
volume, from (x) = to (x) = 262 inches, of the body whose plan and elevation 
contour lines are given by the equations: 

y - N(0.66636)(x)°* 4 -N(0. 28565-5) (x) 2 * 625 

OS 2 05 

W' y = N(0.23829)(x) - N(0.73174-4)(x) ' 

Before setting up the equation in the form of Eq. (164), it is first 
necessary to multiply out the constants and exponents, as follows: 

\ fcv (x)** 1 + Dbl * N(0.666S6 + 0.23829) (x) * 4 + °' 6 

= N(0. 90465) (x) 

n. + n r, x , vl , ,0.6 + 2.625 
k^ k^(x)^l a 2 = N [(0.28565-5) + (0.23829)J (x) 

= N(0.52394-5)(x) 3,225 
k b 1^ (x) n& l * ° b 2 » N [(0.73174-4) + (0.66636)] (x)°* 4 . + 2 *° 5 

= N(0.39810-3)(x) 2 * 45 
k b2 k ag (x) na 2 + ^2 = N [(0.73174-4) + (0.28565-5)] (x) 2 ' 625 + 2 '° 5 

» N(0.01739-8)(x) 4 ' 675 

Substituting these terms in an equation similar to Eq. (165), the 
volume is given by: 



*-> 
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v = 



(TT)N(0.90465)(x) < 



262 



(Tf)N(O. 52594-5) (x) 



4.225 



4.225 



(TT)N(0. 59810-5) (x) 



5.45 



5.45 



262 



J 



(Tr)N(0.01759-8)(x) 



- 262 

^o 

5.675 



5.675 



262 
J 



For (x) ■ 262 inches, this becomes: 



log(V/TT) = 



[(0. 90465) +■ 2 log(262) - log(2)] - [(0.52594-5) + 4.225 log(262) 
- log(4.225)] - [(0.39810-3) + 3.45 log(262) - log(5.45)] + 
[(0.01739-8) + 5.675 log(262) - log(5.675)] 
- 5.44022 -5.1154S - 5.20342 + 4.98727 

V= (275,563-130,445-159,743 + 97,112) 
=259,140 cubic inches 

5.60 - Proportional Change of Volume of a Body of Revolution ;- In many instances, 
after a body of revolution » has been developed from a power curve equation, and 
the volume checked, it will be found that the volume is smaller or larger than 
desired, and a new equation is required. 

Instead of going through all the processes of developing a new equation, 
however, the required equation can be more easily found by employing the proced- 
ure shown below. 

The volume of a body of revolution is equal to the length, (L) , multi- 
plied by the average cross-sectional area, which, being circular, is given by: 

where (A) is the cross-sectional area and (r) the radius. Since, in power 
curve equation form, (r) a is the (7) avg » the volume is: 



I 






WifiW 






' U5. : 
V= (Tr)L(y)J vg ..(166) 

Now, if the letter (f) represents what may be called the "average 
fineness ratio", then: 

f=V(y) avg -c 167 ); 



or: 



L = f fr>avg ....;..........(168) 



Substituting Eq. (168) in Eq. (166): 



>3 



:;. v = (ir)f( 7 ); vg 

or: 

, (y) avg = (v/(Tf)f ) lA . ...(16$) 

If ( v l) represents the volume of the body as first developed, (Vg) the 

desired volume, (yi) and (yo) the respective average ordinates, then: ■' 

avg , *■ avg 

CygW^lW'/^^/CVi/lff!) 175 ...............(170) 

By keeping the "average fineness ratio" of both bodies the same, (f^) 
and (fg) in Eq. (170) cancel out, and (y) maybe substituted for (y) aVff . 
Hence: 

^W^i'-c = ( VV Vs - ••••••• •&«)' 



In a similar manner it can be shown that: 

1/3 v 

V 1 !^ (V 7 !) ••• - ....-.....■.■;:(ire); 



H-,T Eqs. (171) and (172) then show how the maximum ordinates and overall 

it'%. , ' ' '■■■'■■' : ' '■'■ '■•,'' 

lengths of two bodies vary, provided the fineness ratios are the same. The equat- 



ion of the first body can then be corrected by these ratios in the manner out- 






sf*? 
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-lined in Section 2.70. 

For an example, consider again the body of revolution whose volume was 
found in Section 5.40 to be 333,106 cubic inches, and whose equation is; 

7 - N(0.66656)(x)°* 4 - N(0.28565-5)(x) 2 * 625 

The length of this« body is 262 inches, and the maximum ordinate of the 
power curve is 26 inches. (See Section 2.42). Let it now be required to develop 
a new power equation for a similar body with the volume increased to 369,600 
cubic inches. Then, from £q. (171): 

(^W^max = (36960Q/555106) 175 
- N(0. 01505) 

The new maximum ordinate of the power curve will then be: 

^^max = lo e( 26 ) + (0-0 1505 ) 
» 1.41497 + 0.01505 

, - 1.43002 

■ > 

(y ) = 26.917 inches 
<■ max 

From Eq. (172), the correction for length is also the same, and therefore 
the new length is: 

log(Lg) = log(262) + 0.01505 
= 2.41830 +• 0.01505 
» 2.43335 
(Lg) * 271.24 inches 

In Section 2.70 it was shown that, in order to change a power curve equat- 
ion for a change in thickness ratio, or maximum ordinate, it is necessary to 
multiply only the constants of the equation by the ratio of the new to the old 



i 



4Ji 



\ 

■k 



^«E«»l««^|g],-.-i«a^,: ^ 
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ordinates, as the exponents are not affected by changes in thickness. Therefore, 
the equation becomes, corrected for the differences in (y )_-,_. only: 



'max 



.625 



y » N(0.66636 + 0.01505) (x) 0,4 - N [(0.28565-5) + (0.01505)] (x) 2J 

= N(0.68141)(x)°' 4 - N(0. 50070-5) (x) 2 ' 625 ........... (175) 

To correct the values of (x) for the change in length, the multiplier 
is the ratio of the old to the new lengths, (or the reciprocal of the ratio 
used for correcting the constants), raised to the same power as (x). This is 
shown in Section 2.70. In log form, this correction is: 

- N [(1.0000-1) - (0.01505)] 

= n(o.98495-i) ,;■;;■',;■■. '^;v 

Inserting this value in Eq. (173): 

y - n[(0.68141) + (0.98495-l)°- 4 ](x)°- 4 -N [(0.3007-5) + (0.98495-l) 2 * 625 ](x) 2 * 625 

- N [(0.68141) + (0.'39398-0.4) + (0.6 - 0.6)] (x)°* 4 - N [(0.30070-5) + (2.58549 - 
2.625) + (0.375 - 0.375)](x) 2 * 625 

» N(0.67539)(x) 0,4 - N(0.26119-5)(x) 2 * 625 ...(174) 

Eq. (174) is then the final equation for the desired body of revolution - 
with a voluae of 569,600 cubic inches. 

5.70 - Proportional Change of Volume of a Body with Elliptical Cross Section ; - 
The procedure used to correct the contour line equations of a body with an r 
elliptical cross section for a change in volume is substantially the same as 
given in Section 5.60 for a body of revolution. However, in this case, the 
(y) T used represents the average of the (y) max values for the plan view and 

elevation view equations, proportioned in the ratio of the major and minor axes. 



4$ k<$: - :. 
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1118 fr)max values are here designated (y) ma j or and (y) fflinor . 

Thus, if the ratio of the values of (y^joj. to (r) m ij X0T for the two 
contour lines is 1.5 to 1.0, then; 

fr>*ajor/fr>miuor = (l-5)/(l-0) 
and the average (y) max is given by: 

2(y) max = (y) major + <7) minor 

= 1 - 5(y) n,inor + ^roinor 



= 2.5(y) 



minor 



Then: 



(?Ur = 2fr W 2 - 5 



and: 



Kiajor = 2 - 5 H,ax/ 2 
If the ratio of (y) ma j OE . to (y)^ WftT . is called (p), then: 

' (7)minor = 2(7)n,axAP + X > ( 175 ) 

(3T>major = 2 (3 r ) m axAp + X > ( 176 ) 

As an example, consider again the elliptical body whose contour lines 
are given by the equations in Section 5.50: 

(y) , n = N(0.66636)(x)°* 4 - N(0.28565-5)(x) 2 * 625 
(y) minor = N(0.23829)( X )°' 6 - N(0.73174-4)(x) 2 - 05 



.'T 



"'■i-K' 



•"& 



i 



A 






The volume of this body, as found in Section 5.50, was 259,140 cubic 
inches. Its length is 262 inches, and the (y)^ of the najor curve is 26 inches 
and that of the minor curve 20.8 inches. Let us assume that the volume is to be 
changed to 251,000 cubic inches, with the ratio of major to minor axes, and the 
"average fineness ratio", kept unchanged. 

From Eq. (171), the ratio of the average (y) max is: 

1/3 
(7 2 ) /(yi) = (231000/259140) 

4 max max - •. ^ 

= N(0.98356-l) ■: 

The average (y) max of the original body was: 

( T1 ) : = (26.0 + 20.8)/2 ;. , 

VJrj -'max 

■ = 23.4 inches 
and hence the average (y)^ of the second body will be: 

( y > S5 Nflog(23.4) + (0.98358-1)] 



t I 



max 



a N(l. 35258) 
'* » 22.521 inches 



for the length, as given by Eq. (172), 



Using the same correction factor 
the new length becomes: 

L 2 »N[log(262) +(0.98336-1)] 

= N(2.40166) 
= 252.15 inches 

The ratio of (y)^ of the major and minor curves is: 



f 
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. and:- . - • ... 

1/p = 0.8 
Therefore, for the new body, from Eqs. (175) and (176): 

^minor " ^A 1 ' 25 * D(22.521) * 

s 20.018 inches 

(y) maJor * ( 2 )/(°\ 8 *•!) (22.521) 
* 25.023 inches 

Correcting the equations of the new body for the change in ordinates, 
as was done in Section 5.60, we get: 

Tminor = N [(0.23829) + (0.98536-1)] (x)°- S - N [(0.73174-4) 
+ (0.98336-1)] (x) 2 * 05 

« N(0.22165)(x)°' 6 - N(0. 71510-4) (x) 2 * 05 

y . or = H [(0.66636) +' (0.98536-1) ](x) - 4 - N [(0.28565-5) 
+ (0.98356-1)] (x) 2 ' 625 
= N(0.64972)(x)°' 4 - N(0. 26901-5) (x) 2 * 625 

To correct the values of (x) for the change in length, the multiplier 
is again the reciprocal of the correction factor, raised to the exponent of (x), 
which is, in log form: 

• (1.0000-1) - (0.98336-1) 
= 0.01664 

Substituting this, in the above equations for (y) m ^ por and (7) ma ^ 0T t 
there results: 



-■?* 



■'. 
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'<gk$ 



'manor 



= N [(0.22165) + (0.01664)°- 6 l(x) - 6 - N [(0.71510-4) 

+ (0.01664) 2 *° 5 ](x) 2 ' 05 
» N(0.23163)(x) * 6 - N(0. 74923^*) (x) 2 * 05 , 



(177) 



and: 



ymajor = N [(0.64972) + (0.01664)°- 4 ](x) - 4 - N [(0.26901-5) 

\2.625"l, .2.625 
+ (0.01664) J(x) 

s N(0. 65638) (x) * 4 - N (0. 31269-5 )(x) 2 * 62 f ...... 



,(178) 



Eqs. (177) and (178) are then the desired equations to give an elliptical 
body of volume equal to 231,000 cubic inches, -with a length equal to 252.15 
inches, and maximum ordinates, in two planes, of 20.018 and 25.023 inches. 



